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1.  Program  Overview 


1.1.  Background 

In  June  1994,  the  State  Water  Resource  Control  Board  (SWRCB)  contracted  with  the 
Lawrence  Livermore  National  Laboratory/  University  of  California  (LLNL/UC)  Leaking 
Underground  Fuel  Tank  (LUFT)  Team  to  study  the  cleanup  of  LUFTs  in  California.  The  study 
consisted  of  data  collection  and  analysis  from  LUFT  cases  and  a  review  of  other  studies  on  LUFT 
cleanups.  Two  final  reports  were  submitted  to  the  SWRCB  in  October  and  November  1995. 
These  reports  were  entitled:  Recommendations  To  Improve  the  Cleanup  Process  for  California’s 
Leaking  Underground  Fuel  Tanks  (LUFTs),  and  California  Leaking  Underground  Fuel  Tank 
(LUFT)  Historical  Case  Analysis. 

Prior  to  submittal,  both  LLNL/UC  reports  were  subjected  to  extensive  internal  peer  review 
within  the  LLNL,  as  well  as  within  each  of  the  four  UC  campuses  involved:  UC  Berkeley, 
UC  Santa  Barbara,  UC  Davis,  and  UC  Los  Angeles.  The  U.S.  Environmental  Protection  Agency 
(EPA)  has  reviewed  both  LLNL/UC  reports  and  issued  a  fact  sheet  supporting  the  findings  and 
recommendations.  Both  reports  have  received  national  attention  and  have  been  extensively 
reviewed  by  many  interested  parties. 

Data  were  collected  primarily  from  the  alluvial  geologic  settings  typical  of  the  San  Francisco 
Bay  Area,  Los  Angeles  Basin,  and  the  Central  Valley.  These  data  represent  California’s  highly 
populated  areas  where  most  gas  stations  are  found  and  most  LUFT  releases  occur.  Study  results 
can  be  applied  to  similar  settings  throughout  California.  The  length  of  dissolved  benzene  plumes 
in  groundwater  at  LUFT  releases  were  evaluated  over  time  to  determine  how  plumes  behave.  Over 
1 ,200  LUFT  cases  were  evaluated.  Benzene  was  analyzed  because  it  is  the  human  carcinogen  of 
greatest  concern  in  fuel.  It  is  relatively  soluble  in  water,  and  cleanup  standards  are  generally  tied  to 
benzene  concentrations. 

Ninety  percent  of  the  dissolved  benzene  plumes  were  less  than  280  ft  in  length.  Most  of  these 
plumes  were  either  stable  or  shrinking  in  length.  Seventy  percent  of  the  plumes  in  the  study  sites 
were  found  in  shallow  groundwater,  less  than  25  ft  below  the  ground  surface. 

1.1.1.  The  Use  of  Natural  Attenuation 

The  study  concluded  that  with  rare  exceptions,  petroleum  fuel  releases  will  naturally  degrade 
(passive  bioremediation)  in  California’s  subsurface  conditions.  Removing  the  source  of  the  release 
will  speed  the  cleanup  time.  Source  removal  includes  removing  leaking  tanks  and  lines,  and 
removing  free  product  and  petroleum  fuel  saturated  soil,  as  much  as  economically  and  technically 
feasible.  The  use  of  passive  bioremediation  still  requires  site  characterization,  an  assessment  of 
potential  risks,  and  an  evaluation  of  the  applicability  of  passive  bioremediation.  Effective  site 
cleanup  can  be  achieved  when  passive  bioremediation  is  appropriately  used. 

Passive  bioremediation  may  be  evaluated  through  both  primary  or  secondary  field  evidence. 
Primary'  evidence  includes  quantitative  evaluation  of  plume  stability  or  plume  shrinkage  based  upon 
trends  in  historical  groundwater  quality  data.  Secondary  evidence  includes  indirect  indicators  of 
bioattenuation,  such  as  variations  in  key  geochemical  parameters  (dissolved  oxygen,  nitrate, 
sulfate,  iron,  manganese,  methane,  alkalinity/carbon  dioxide,  Eh,  pH)  between  measurements  in 
hydrocarbon-impacted  areas  and  background.  Although  primary  evidence  generally  provides  the 
strongest  arguments  to  support  natural  attenuation  at  a  given  site,  such  evidence  is  commonly  not 
available  because  adequate  historical  groundwater  monitoring  may  not  exist.  Consequently,  means 
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for  assessing  the  role  of  natural  attenuation  in  controlling  risk  by  secondary  evidence  must  be  more 
fully  explored. 

1.2.1.  Risk-Based  Corrective  Action 

The  study  also  concluded  that  risk-based  corrective  action  (RBCA)  provides  a  framework  to 
link  cleanup  decisions  to  risk.  For  a  risk  to  exist,  there  must  be  a  source  of  a  hazard,  a  receptor, 
ar)da  path  way  that  connects  the  two.  All  three  factors  must  be  addressed  to  determine  whether  a 
LUFT  release  poses  a  risk  to  human  health,  safety,  or  the  environment.  If  the  source,  pathway,  or 
receptor  are  at  all  times  absent,  there  is,  by  definition,  no  risk. 

Risk  characterization  is  defined  as  an  information  synthesis  and  summary  about  a  potentially 
hazardous  situation  that  addresses  the  needs  and  interests  of  decision  makers  and  of  interested  and 
affected  parties.  Risk  characterization  is  a  prelude  to  cleanup  decision  making  and  depends  on  an 
iterative,  analytic,  and  deliberative  process.  This  process  attempts  to  gather  all  relevant  data  so  the 
decision  makers  may  then  choose  the  best  risk-management  approach. 

If  the  risk  of  affecting  receptors  (humans  or  ecosystems)  is  low,  then  the  following  risk 
management  strategy  is  appropriate  and  cost  effective:  1)  remove  and  upgrade  leaking 
underground  tanks  and  lines,  2)  perform  source  removal  as  much  as  economically  and  technically 
feasible,  and  3)  use  passive  bioremediation  at  low-risk  sites,  supported  by  monitoring  as 
appropriate. 

1.2.  LUFT  Demonstration  Program 

One  of  the  important  recommendations  of  this  study  was  to  identify  a  series  of  LUFT 
demonstration  sites  and  to  form  a  panel  of  experts  made  up  of  scientific  professionals  from 
universities,  private  industry,  and  Federal  and  State  regulatory  agencies.  This  panel  would  provide 
professional  interpretations  and  recommendations  regarding  LUFT  evaluations  and  closures  at 
demonstration  sites. 

As  a  result  of  this  recommendation,  nine  Department  of  Defense  (DoD)  sites  were  selected. 
Site  selection  was  coordinated  through  the  California  Military  Environmental  Coordination 
Committee  (CMECC)  Water  Process  Action  Team  (PAT).  Sites  were  selected  to  represent  each 
branch  of  the  military  services  with  bases  in  California,  as  well  as  a  number  of  Regional  Water 
Quality  Control  Boards  (RWQCB)  and  the  diverse  hydrogeologic  settings  in  California  where  fuel 
hydrocarbon  contaminant  (FHC)  cleanup  problems  occur.  The  sites  selected  and  their 
corresponding  RWQCB  region  are: 

•  Army  Presidio  at  San  Francisco,  San  Francisco  RWQCB. 

•  Barstow  Marine  Corps  Logistic  Center,  Lahonton  RWQCB. 

•  Castle  Air  Force  Base,  Central  Valley  RWQCB. 

•  China  Lake  Naval  Weapons  Center,  Lahonton  RWQCB. 

•  El  Toro  Marine  Corps  Air  Station,  San  Diego  RWQCB. 

•  George  Air  Force  Base,  Lahonton  RWQCB. 

•  Port  Hueneme  Naval  Construction  Battalion  Center,  Los  Angeles  RWQCB. 

•  Travis  Air  Force  Base,  San  Francisco  RWQCB. 

•  Vandenberg  Air  Force  Base,  Central  Coast  RWQCB. 
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The  Expert  Committee  (EC)  selected  to  evaluate  the  selected  demonstration  sites  are: 

•  Mr.  David  W.  Rice,  LLNL,  Environmental  Scientist;  Project  Director  SWRCB  LUFT  Re- 
evaluation  Project;  LLNL/UC  LUFT  Team  member;  DoD  FHC  Demonstration  Program 
Coordinator. 

•  Dr.  Walt  McNab,  LLNL,  Hydrogeochemist,  with  expertise  in  the  evaluation  of  passive 
bioremediation  processes. 

•  Dr.  William  E.  Kastenberg,  UC,  Berkeley,  Professor  and  Chairman,  Department  of 
Nuclear  Engineering;  member  of  LLNL/UC  LUFT  Team,  with  expertise  in  environmental 
decision  making  and  decision  analysis  processes. 

•  Dr.  Lome  Everett,  UC,  Santa  Barbara,  Hydrogeologist;  head  of  the  Vadose  Zone  Research 
Laboratory  and  member  of  LLNL/UC  LUFT  Team,  with  expertise  in  vadose  zone  PHC 
transport  mechanisms  and  passive  bioremediation  processes. 

•  Dr.  Stephen  Cullen,  UC,  Santa  Barbara,  Hydrogeologist;  member  of  LLNL/UC  LUFT 
Team  with  expertise  in  vadose  zone  PHC  transport  mechanisms  and  passive  bioremediation 
processes. 

•  Dr.  Paul  Johnson,  Arizona  State  University,  Chemical  Engineer;  primary  author  of 
American  Society  for  Testing  and  Materials  (ASTM)  RBCA  guidance,  with  expertise  in 
chemical  fate  and  transport. 

•  Dr.  Michael  Kavanaugh,  Chairman,  National  Research  Council  Alternatives  for 
Groundwater  Cleanup  Committee;  Environmental  Engineer,  Environ  Corporation,  with 
expertise  in  evaluation  of  groundwater  remediation  alternatives  and  environmental  decision 
making  processes. 

•  Mr.  Matthew  Small,  U.S.  EPA  Region  IX,  Hydrogeologist;  Co-Chairman  of  U.S.  EPA 
Remediation  by  Natural  Attenuation  Committee,  with  expertise  in  risk-based  corrective 
action  and  passive  bioremediation. 

The  demonstration  program  process  can  be  summarized  in  the  following  nine  steps: 

Step  1:  Site  scoping  meeting  with  site  staff,  regulators,  and  EC  staff  representatives.  Develop 
and  discuss  site  conceptual  model.  Identify  and  discuss  pathways  and  receptors  of  concern. 

Step  2:  RBCA  training  for  DoD  Petroleum  Hydrocarbon  Cleanup  Demonstration  Program 
(PHCDP)  participants. 

Step  3:  Site  staff  and  contractors  prepare  the  data  package.  EC  staff  reviews  available  data  and 
identifies  data  gaps  needed  to  apply  RBCA  approach. 

Step  4:  EC  visits  site  and  receives  briefing,  on  site  characterization,  conceptual  model,  and 
pathways  and  receptors  of  concern.  Site  tour  is  included  in  this  briefing.  Following  EC’s  visit,  a 
site  characterization  report  is  prepared  by  the  EC  containing  recommendations  for  further  data 
collection,  if  needed. 

Step  5:  EC  staff  applies  RBCA  approach  to  the  site  using  best  available  data. 

Step  6:  EC  staff  evaluates  the  natural  attenuation  potential  for  the  site  using  best  available  data. 
An  estimate  of  the  time  to  clean  up  and  the  uncertainty  associated  with  this  estimate  will  be  made. 
Sampling  and  monitoring  procedures  to  support  intrinsic  bioremediation  for  the  site  will  be 
identified. 
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Step  7:  Based  on  the  concept  of  applied  source,  pathways,  and  receptors  as  to  potential 
hazards,  site  specific  findings  regarding  natural  attenuation  potential,  and  discussion  with 
regulators,  the  EC  shall  provide  its  recommendations  for  an  appropriate  risk-management  strategy 
at  the  site  and  the  set  of  actions  needed  to  achieve  site  closure.  The  EC  will  present  its 
recommendations  at  an  appropriate  forum. 

Step  8:  The  EC  will  provide  aDoD  PHCDP  overall  evaluation  comparing  the  effectiveness  of 
CalRBCA  at  each  site  in  the  program.  An  estimation  of  the  cost  savings  using  the  CalRBCA 
protocol  will  be  compared  to  baseline  approaches.  An  estimation  of  the  value  of  the  remediated 
water  will  be  made. 

Step  9:  The  EC  Staff  will  produce  a  DoD  Risk  Execution  Strategy  for  Clean-Up  of  the 
Environment  (RESCUE)  implementation  guide  and  accompanying  procedures  manual  (Phase  I, 
Petroleum)  that  can  be  used  in  California  and  in  other  states  by  military  bases. 


2.  Site  Overview 


2.1.  Background 

George  Air  Force  Base  (GAFB)  is  located  in  Victorville,  California.  The  base  was  active  until 
1992  when  it  was  closed  as  part  of  the  base  realignments  and  closures  approved  by  the  Secretary 
of  Defense  and  Congress.  It  is  anticipated  that  GAFB  will  eventually  be  used  as  a  regional 
international  airport.  The  base  facilities  are  intended  to  house  a  variety  of  small  businesses  and 
light  industry. 

Petroleum  hydrocarbons  (originally  JP-4  aviation  fuel)  are  in  the  subsurface  environment  at 
GAFB  as  a  result  of  leaks  from  underground  tanks  and  piping  systems.  A  significant  volume  of 
separate  phase  product,  estimated  to  be  on  the  order  of  350,000  gal,  persists  on  top  of  the  water 
table,  covering  an  area  of  approximately  32  acres,  with  product  thickness  ranging  from  a  thin 
sheen  to  several  feet.  A  plume  of  dissolved  FHCs  exists  in  the  underlying  groundwater,  extending 
several  hundred  feet  downgradient  of  the  free  product  lens.  Isolated  pockets  of  low  concentrations 
of  dissolved  halogenated  hydrocarbons  (trichloroethylene  [TCE],  tetrachloroethylene  [PCE], 
dichloroethane  [1,1-DCA])  have  also  been  noted  in  site  groundwater. 

GAFB  Operable  Unit  #2  (OU#2)  was  selected  as  a  site  for  the  Demonstration  Program  because 
the  base  has  collected  a  considerable  amount  of  data  and  has  established  technical  protocols  that 
would  be  valuable  in  evaluating  the  recommendations  proposed  by  the  LLNL/UC  report. 
Extensive  remedial  investigations  and  risk  assessments  have  been  performed  according  to  EPA’s 
Risk  Assessment  Guidance  for  Supetfund  -  Human  Health  Evaluation  Manual  and  were  available 
for  review.  The  petroleum  hydrocarbon  compounds  benzene,  toluene,  ethylbenzene,  and  xylenes 
(BTEX)  associated  with  OU#2  were  identified  as  chemicals  of  concern. 

2.2.  Site  Conceptual  Model 

The  site  is  underlain  by  alluvial  sediments  consisting  of  unconsolidated  sand,  silt,  clay,  and 
gravel  in  an  arid  climatic  setting.  A  perched  unconfined  aquifer  extends  from  a  depth  of  about  125 
to  196  ft  and  is  perched  on  a  20-ft-thick  aquiclude  of  clay  and  clayey  silt.  Below  the  perched 
aquifer  is  a  second  vadose  zone  (unsaturated  zone)  extending  to  a  deeper  aquifer  at  a  depth  of 
268  ft.  The  aquiclude  is  continuous  from  exposures  southeast  of  the  base  to  north  of  OU#2. 

Hydraulic  conductivity  values  have  been  calculated  from  a  number  of  pumping  and  aquifer 
tests  (although  analyses  associated  with  the  two  methods  have  produced  estimated  values  that  are 
not  entirely  self-consistent).  Taken  with  interpolated  hydraulic  head  values  and  porosity  estimates, 
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these  data  suggest  groundwater  movement  in  a  north-northeasterly  direction  at  a  rate  of  tens  to 
hundreds  of  feet  per  year.  Little  recharge/evapotranspiration  is  expected  in  association  with  the 
unconfined  aquifer. 

The  spatial  extent  of  the  dissolved  FHC  constituents/BTEX  in  groundwater  appears  reasonably 
well-characterized  over  individual  sampling  rounds,  although  long-term  historical  trends  in 
individual  wells  are  not  available.  Concentrations  of  natural  attenuation  geochemical  indicator 
species  (e.g.,  dissolved  oxygen,  nitrate,  sulfate,  iron,  methane,  redox  potential)  have  also  been 
collected  from  16  monitoring  wells  as  part  of  a  natural  attenuation  monitoring  study.  Clear 
evidence  has  been  presented  indicating  natural  bioattenuation  processes  are  occurring  at  GAFB. 
These  include  apparent  shrinkage  of  BTEX  constituent  plumes  with  time,  loss  of  BTEX 
components  in  comparison  to  recalcitrant  dissolved  FHC  constituents  (trimethylbenzene  isomers), 
and  changes  in  geochemical  indicator  species  which  vary  in  space  in  a  manner  consistent  with 
expected  bioattenuation. 

A  review  of  the  adequacy  of  the  site  conceptual  model  was  presented  to  GAFB  in  a  letter  dated 
July  30,  1996  (Appendix  A). 


3.  Risk  Analyses 


3.1.  Sources 


3.1.1.  Primary  Source(s) 

Petroleum  hydrocarbons  (originally  JP-4  aviation  fuel)  are  in  the  subsurface  environment  at 
GAFB  as  a  result  of  leaks  from  underground  tanks  and  piping  systems.  Vadose  zone 
contamination  encountered  during  drilling  and  soil  analyses  indicate  that  there  were  two  primary 
sources  for  the  free  product  plume:  Fuel  Pit  1  and  Fuel  Pit  3/Main  Pump  Station  area  located  in  the 
flight-line  aircraft  parking  apron.  The  original  liquid  fuel-supply  system,  constructed  of  3-  to  6-in. 
aluminum  pipes,  leaked  an  unknown  volume  during  the  early  1970s.  Four  50,000-gal 
underground  storage  tanks  (USTs)  were  taken  out  of  service  in  1990  after  failing  two  consecutive 
leak  tests.  The  rates  and  quantities  of  leakage  are  unknown.  The  USTs  were  removed  along  with 
contaminated  soils  and  disposed  of. 

3.1.2.  Secondary  Source(s) 

A  significant  volume  of  separate-phase  product,  estimated  on  the  order  of  350,000  gal,  persists 
on  top  of  the  perched  aquifer  water  table.  This  free  product  covers  an  area  of  approximately 
32  acres,  with  product  thickness  ranging  from  a  thin  sheen  to  several  feet.  The  free  product  plume 
appears  to  flow  to  the  northwest  from  release  areas  under  the  flight-line  aircraft  parking  apron. 
Movement  of  the  free  product  is  slowed  and  possibly  stopped  by  the  fine-grained  deposits  within 
buried  stream  channels.  No  contamination  has  been  detected  in  groundwater  samples  taken  from 
the  deep  aquifer  beneath  OU#2. 

3.1.3.  Engineered  Source  Removal 

Current  active  source  remediation  techniques  in  use,  or  under  consideration,  at  GAFB  include 
free  product  removal  and  capillary  fringe  bioventing.  Free  product  remediation  consists  of  mobile 
and  permanent  free  product  recovery  systems,  which  remove  approximately  160  gal  of  free 
product  per  day,  with  a  cumulative  total  of  approximately  32,000  gal  to  date.  The  bioventing 
program  has  to  date  consisted  largely  of  a  treatability  test. 
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Experience  suggests  that  free  product  removal  tends  to  be  less  successful  the  longer  the  time 
since  release.  Based  on  current  free  product  estimates,  the  current  mass  removed  represents 
approximately  10  percent  of  the  estimated  free  product  source  mass. 

For  new  releases,  removal  volumes  between  70  percent  and  90  percent  are  possible,  depending 
on  the  geological  conditions  at  the  site.  This  percentage  quickly  declines  as  the  non-aqueous  phase 
liquid  (NAPL)  spreads,  becoming  sorbed,  trapped  by  capillary  forces,  and  isolated  as  distinct 
ganglia.  Removals  of  up  to  30  percent  are  much  more  typical. 

3.1.4.  Passive  Source  Removal  By  Natural  Dissipation  Processes 

To  address  the  time  required  for  the  free  product  lens  to  dissolve  (coupled  with  bioattenuation 
in  the  dissolved  phase),  IT  Corporation  developed  a  site  plume  model  using  the  BIOPLUME  II 
simulator  to  predict  the  behavior  of  the  dissolved  hydrocarbon  plume,  as  influenced  by 
bioattenuation.  This  plume  model  accounted  for  contributions  of  fuel  hydrocarbons  to  the 
dissolved  phase  by  dissolution  of  the  separate-phase  source  term  over  time.  The  results  of  the 
BIOPLUME  II  simulations  suggested  that  the  separate-phase  source  would  disappear  after 
approximately  20  years,  with  the  dissolved  plume  persisting  for  an  additional  20  years  until 
bioattenuation  reduced  concentrations  below  Maximum  Contaminant  Levels  (MCLs). 

In  our  opinion,  the  estimated  20-year  dissolution  period  calculated  in  the  original 
BIOPLUME  II  study  cannot  be  substantiated  because  of  problematic  assumptions  in  the  source 
diffusion  model.  Our  calculations,  using  two  alternative  modeling  approaches  that  take  into 
account  solubility  limitations  of  hydrocarbon  compounds  in  water,  suggest  much  longer  separate- 
phase  source  dissolution  times,  on  the  order  of  50  to  100  years.  However,  we  hypothesize  that 
volatile  losses  of  hydrocarbon  components  into  the  vadose  zone,  a  process  not  considered  in  the 
initial  site  BIOPLUME  II  modeling,  may  represent  a  far  more  significant  mechanism  for  source 
removal,  potentially  resulting  in  the  depletion  of  the  separate-phase  source  in  approximately  one- 
half  the  estimated  dissolution  time  (i.e.,  a  “best-case”  estimate  on  the  order  of  30  years). 
Nevertheless,  there  is  a  high  degree  of  uncertainty  associated  with  this  volatilization  analysis 
because  of  the  large  number  of  poorly  understood  factors  that  influence  transport  of  volatiles 
through  the  vadose  zone  (e.g.,  spatial  variability  in  soil  moisture  content,  lithology,  transient 
nature  of  the  separate-phase  source,  barometric  pumping,  etc.).  As  a  result,  the  separate-phase 
source  depletion  estimates,  combined  with  the  BIOPLUME  II  dissolved  phase  bioremediation 
estimates,  cannot  produce  a  precise  calculation  for  time  to  achieve  cleanup.  In  our  opinion, 
remediation  of  the  dissolved  hydrocarbon  plume  at  GAFB  is  probably  not  likely  to  occur  in  a  time 
frame  of  less  than  50  years. 

Calculations  pertaining  to  dissolution  and  volatilization  of  the  free  product  lens  are  presented  in 
Appendix  B. 

322.  Pathways  and  Fate 

To  pose  a  human  health  or  ecological  risk,  the  source  of  contaminants  must  be  linked  to 
receptors  via  pathways.  In  the  case  of  groundwater,  advective,  dispersive,  or  diffusive  transport 
of  dissolved  contaminants  to  a  water-supply  well  would  constitute  such  a  pathway.  However, 
bioattenuation  processes  tend  to  limit  the  migration  of  dissolved  hydrocarbons;  the  input  of 
dissolved  contaminants  from  residual  sources  is  balanced  by  bioattenuation  integrated  over  the 
extent  of  the  plume.  This  mass  balance  constraint  is  the  likely  explanation  for  the  limited 
dissolved-phase  plume  lengths  associated  with  majority  of  LUFT  sites  (Rice  et  al„  1995).  Thus, 
the  problem  of  assessing  the  likelihood  for  contaminant  migration  in  groundwater  is  linked  directly 
to  bioattenuation. 

Data  collected  from  the  site  provide  several  lines  of  evidence  that  support  biodegradation 
reactions  of  fuel  hydrocarbons.  Changes  in  concentration  of  various  BTEX  constituents  with  time 
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are  the  least  reliable  among  these  lines  of  evidence  because  of  the  relatively  short  period  during 
which  the  site  has  been  monitored  relative  to  anticipated  rates  of  changes  in  concentration,  and  the 
presumed  persistence  of  the  separate-phase  source  term.  Changes  in  the  concentration  ratios  of 
BTEX  components  and  recalcitrant  tracers  such  as  trimethylbenzene  do  provide  additional  evidence 
supporting  bioattenuation.  However,  in  the  absence  of  historical  data  directly  demonstrating  plume 
attenuation,  quantitative  analysis  of  plume  behavior  under  the  influence  of  bioattenuation  may 
instead  incorporate  predictive  modeling  or  analysis  of  secondary  evidence. 

3.2.1.  Groundwater  Plume  Stability  Analysis 

At  the  GAFB  site,  great  uncertainty  exists  concerning  the  length  of  time  required  for  depletion 
of  the  separate-phase  source.  Thus,  a  conservative  model  assuming  constant  source  is  warranted 
for  forecasting  plume  behavior.  Under  such  conditions,  the  plume  concentration  would  be 
expected  to  stabilize  in  a  steady-state  condition.  Using  a  simple  analytical  solution  to  the 
advective-dispersive  transport  problem  (Domenico,  1987),  probabilistic  modeling  was  conducted 
to  forecast  plume  lengths  based  upon  site-specific  probability  distributions  of  governing  parameters 
(e.g.,  hydraulic  conductivity,  first-order  hydrocarbon  decay  rate,  length  of  time  source  has  been 
active).  Thus,  although  sufficient  historical  monitoring  data  are  lacking,  a  simple  model  of 
expected  plume  behavior  can  be  constructed  to  quantify  plume  length  as  a  function  of  time. 

Benzene  plume  lengths  at  GAFB,  based  on  this  type  of  modeling,  generally  match  the 
observed  plume  length.  Projection  of  plume  behavior  into  the  future  (i.e.,  to  the  year  2022),  under 
the  assumption  of  a  constant  source  (a  conservative  assumption)  and  constant  hydrogeologic 
parameters,  suggests  that  the  probability  that  the  plume  will  advance  over  significantly  greater 
distances  is  very  low.  Thus,  taking  into  account  die  length  of  time  the  contaminant  source  has 
been  in  existence,  it  is  likely  that  the  plume  is  approaching,  or  has  already  achieved,  a  steady-state 
condition.  Mass  balance  calculations  are  also  consistent  with  this  conclusion. 

The  calculations  used  in  the  modeling  are  presented  in  Appendix  C. 

3.2.2.  Evaluation  of  Geochemical  Indicators  of  Bioattenuation 

Geochemical  indicators  of  biodegradation  provide  evidence  that  is  less  dependent  on  the  length 
of  the  monitoring  period.  Concentrations  of  natural  attenuation  geochemical  indicator  species 
(e.g.,  dissolved  oxygen,  nitrate,  sulfate,  iron,  methane,  redox  potential)  have  been  collected  from 
16  monitoring  wells  at  GAFB  as  part  of  a  natural  attenuation  monitoring  study.  We  have  evaluated 
the  observed  variations  in  hydrocarbon  concentrations  and  key  geochemical  indicators  of 
biodegradation  across  the  site  using  a  reactive  transport  model  screening  tool  coupled  with  the 
analytical  solute  transport  model  of  Domenico  (1987).  The  purpose  of  the  modeling  was  to  verify, 
in  a  semi-quantitative  manner,  that  the  observed  variations  in  site  geochemistry  are  indeed 
consistent  with  those  expected  as  a  result  of  biodegradation  processes. 

The  modeling  analyses  is  explained  in  Appendix  C.  The  results  indicate  that  biodegradation 
reactions  fully  explain  the  observed  changes  in  groundwater  geochemistry,  supporting  the 
probability  that  natural  attenuation  processes  are  actively  remediating  the  hydrocarbon  plume. 

3.2.3.  Potential  Human  Exposure  Pathways 

Ingestion.  Due  to  the  depth  to  the  free-product  contaminated  portion  of  the  subsurface  (i.e., 
the  water  table),  on  the  order  of  100  ft,  and  the  90-percent  pavement  cover  over  the  site,  the 
ingestion  pathway  is  not  considered  likely.  Part  of  assessing  the  pathway  involves  the  idea  of 
awareness  of  a  hazard.  When  there  is  awareness,  hazards  are  typically  avoided,  thus,  removing 
the  potential  impact  of  the  hazard  and  removing  the  risk.  A  strong  indicator  of  hazard  is  the  low 
taste  and  odor  thresholds  associated  with  FHC  products.  At  levels  where  there  would  exist  a  large 
risk  (10"4  excess  cancer  risk)  the  odors  and  tastes  associated  would  cause  a  voluntary  avoidance  of 
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that  water  for  drinking  and  showering,  although  workers  may  not  show  the  same  avoidance  for 
inhalation  of  fugitive  vapors.  Therefore,  the  levels  associated  with  the  separate  phase  product  at 
GAFB  produce  an  upper  bound  to  the  possible  risk  associated  with  the  site.  Due  to  this,  ingestion 
or  inhalation  of  water  from  showering  are  improbable  pathways,  and  therefore  the  risk  associated 
with  direct  ingestion  or  showering  with  these  waters  is  minimal. 

Inhalation.  Vapor  exposure  pathways  associated  with  the  free  product  lens  are  unlikely  to  pose 
any  substantial  risks  through  inhalation,  due  partly  to  the  large  asphalt/concrete  cap,  and  due  partly 
to  attenuation  in  the  vadose  zone.  Risk  calculations  pertaining  to  the  vapor  exposure  pathway  are 
presented  in  Appendix  D. 

3.3.  Receptors 


3.3.1.  Present  and  Anticipated  Future  Human  Health  Receptors 

Population  centers  are  clustered  south  of  GAFB.  Population  counts  within  five  miles  of  the 
plume  center,  including  6,815  individuals  that  were  part  of  the  GAFB  census  tract,  are  shown  on 
Table  1 .  It  is  assumed  that  this  particular  population  will  diminish  as  GAFB  is  closed,  but  the 
general  population  in  the  area  may  increase  based  on  transformation  of  GAFB  to  the  Southern 
California  International  Airport  (SCIA).  Currently  this  is  only  in  the  planning  stages.  The 
majority  of  the  population  is  located  south  of  GAFB;  to  the  north-east  a  population  of  1,169 
individuals  are  spread  out  over  a  large  area.  These  individuals  are  likely  to  use  private  well  water. 


Table  1.  Population  distribution  around  GAFB. 


Radius  from 

FHC  plume  center 
(miles) 

Population 

1 

6,815 

2 

7,829 

3 

12,890 

4 

15,129 

5 

20,855 

Population  in  Victorville  has  experienced  an  approximate  9  percent/year  exponential  increase 
since  1980.  This  growth  has  decreased  substantially,  and  was  no  more  than  3.5  percent/year  over 
the  last  year.  Figure  1  shows  population  growth,  based  on  various  growth  rates  for  the  Victorville 
area.  Growth  is  greatest  around  the  city  itself,  and  becomes  less  to  the  north.  Figure  2  shows  the 
expected  growth  rate  for  the  total  area.  A  moderate  6  percent  growth  rate  will  increase  water  usage 
approximately  4.5  times.  The  Department  of  Water  Resources  expects  an  overall  increase  of  urban 
growth  demand  of  3  percent/year,  which  may  be  inappropriately  low  considering  the  recent 
growth  of  the  Victorville  and  Antelope  Valley  areas  and  the  possible  conversion  of  GAFB  into  a 
regional  airport.  Groundwater  overdraft  for  the  entire  South  Lahontan  Region  is  approximately 
70,000  acre-ft/year,  which  is  expected  to  remain  constant  due  to  a  reduction  in  agricultural  usage 
and  importing  water  from  the  California  Water  System. 

Victor- Valley  Water  (WW)  expects  to  supply  the  SCIA  with  water  from  a  groundwater  source 
that  will  be  piped  in  from  the  Victorville  area  (Don  Songer,  Personal  Communication,  1996). 
WW  currently  does  not  plan  to  use  the  groundwater  on  the  AFB  site  itself.  The  California 
Department  of  Health  Services  database  shows  supply  wells  within  a  two-mile  radius  of  GAFB, 
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supplying  the  population  to  the  south  of  the  base  and  the  City  of  Adalento.  The  Department  of 
Water  Resources  (DWR)  database  states  that  the  production  capacities  of  these  wells  range  from 
100  to  1,200  gallons  per  minute  (gpm),  and  are  screened  in  the  Regional  (lower)  aquifer  at  depths 
ranging  from  500  to  610  ft  below  ground  surface.  These  wells  are  all  upgradient  of  the  base,  and 
would  not  be  affected  by  any  groundwater  impact.  Downgradient  of  the  site,  Victor  Valley,  Waste 
Water  Reclamation  Authority,  uses  two  wells  screened  in  the  upper  aquifer,  which  have  pumping 
capacities  of  500  gpm.  These  wells  are  used  for  non-potable  industrial  applications  only  (JMMCE, 

Two  miles  north  east,  there  are  about  50  to  60  private  wells,  situated  mostly  along  Highway  66 
(Jeff  Ricciardi,  personal  communication,  1996).  These  wells  draw  their  water  directly  from  the 
Mojave  River  alluvial  bed.  There  are  a  few  wells  between  GAFB  and  the  Mojave  River  bed  area, 
but  the  zone  of  influence  for  such  small  wells  will  not  affect  the  FHC  plume  at  the  base.  There  are 
no  private  water-supply  wells  within  one-half  mile  of  the  release. 

The  potential  for  human  receptors  to  ingest  groundwater  impacted  by  the  FHC  plume  is 
negligible,  considering  current  population  patterns,  proposed  water  plans  for  the  area,  and  the 
apparent  stability  of  the  GAFB  hydrocarbon  plume.  Therefore,  groundwater  is  not  a  health  risk 
pathway. 

No  residences  or  buildings  would  be  affected  by  petroleum  vapors  rising  through  the  soil 
column  due  to  the  FHC  plume  under  the  airfield. 

3.3.2.  Ecological  Risk  Receptors 

No  ecological  risk  receptors  have  been  identified. 


4.  Summary  and  Recommendations 

For  a  dissolved  hydrocarbon  plume  to  pose  a  risk  to  human  health  or  the  environment,  a 
source,  pathway,  and  receptor  must  all  simultaneously  co-exist.  The  absence  of  any  one  of  the 
three  implies  that  there  is  no  associated  risk.  At  GAFB,  each  of  these  components  has  been 
evaluated  independently: 

Sources.  An  estimated  350,000  gal  of  free  product  remain  in  the  subsurface.  It  is  unlikely  that 
current  free  product  removal  efforts  will  be  able  to  remove  the  majority  of  the  contaminants. 
Natural  dissipation  processes,  including  dissolution  and  volatilization,  will  probably  require  at  least 
30  years  to  remove  a  majority  of  the  volatile,  soluble  components  (BTEX). 

Pathways.  In  our  judgment,  historical  groundwater  monitoring  at  GAFB  is  insufficient  to 
provide  clear  primary  evidence  of  plume  stability  or  decline.  Nevertheless,  geochemical  indicator 
data  and  recalcitrant  tracer  analyses  do  provide  secondary  evidence  strongly  supporting  some 
degree  of  bioattenuation.  Contaminant  transport  modeling  suggests  a  high  probability  that  the 
dissolved  plume  has  at  least  achieved  a  steady-state  condition. 

Receptors.  Accounting  for  current  and  projected  future  local  population  and  land  use  issues,  it 
appears  that  there  are  no  receptors  at  risk  of  detrimental  impact  by  the  hydrocarbons  in  the 
subsurface.  Furthermore,  because  the  hydrocarbon  plume  remains  on  Air  Force  property, 
institutional  controls  can  be  implemented  to  assure  against  hydrocarbon  exposure  risks  associated 
with  unanticipated  land  use  or  groundwater  use  activities. 

Taste  and  odor  thresholds  are  a  potent  indicator  of  hazard.  The  presence  of  low  taste  and  odor 
thresholds  associated  with  the  relatively  high  concentrations  seen  in  the  separate  phase  FHC 
regions,  is  a  potent  warning  that  will  cause  potential  receptors  to  avoid  drinking  or  inhaling  during 
showering  these  waters,  and  act  as  an  upper-bound  to  the  potential  risks  associated  with  this  site. 
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In  summary ,  despite  the  presence  of  a  continuing  source,  no  at-risk  receptors  can  be  identified. 
Furthermore,  the  existence  of  migration  pathways  to  carry  the  dissolved  plume  off  the  GAFB  site 
is  highly  questionable  because  the  plume  is  essentially  in  a  steady-state  condition  (or  soon  will  be). 
As  a  result,  the  risk  associated  with  the  hydrocarbon  plume  at  GAFB  is  minimal.  Given  this  low 
risk,  we  believe  that  remediation  by  natural  attenuation,  as  advocated  by  GAFB  and  IT  Corp.,  will 
protect  human  health  and  the  environment  and  is  thus  an  acceptable  means  of  remediation  at  the 
site.  As  a  safeguard  measure,  we  strongly  recommend  that  at  least  two  wells  located  near  the 
downgradient  edge  of  the  dissolved  plume  be  chosen  as  sentry  wells,  and  be  monitored  for  BTEX 
components  annually.  It  is  unclear  how  significant  the  impact  of  continued  free  product  removal 
efforts  may  be  on  eventual  cleanup  of  the  hydrocarbon  contamination.  Given  the  large  quantity  of 
free  product,  it  is  probable  that  die  free  product  recovery  impact  on  the  dissolved  hydrocarbon 
plume  will  be  small  in  comparison  to  natural  attenuation  mechanisms  over  time. 
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Figure  1.  Population  growth  in  Victorville  assuming  exponential  growth. 
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Appendix  A 

George  Air  Force  Base  Site  Assessment 
Evaluation  Letter 


Lawrence  Livermore  National  Laboratory 


July  30, 1996 


Bob  Sommer 
OL-C/AFBCA,  Bldg  321 
13436  Sabre  St. 

Victorville,  CA  92394-5000 

Dear  Bob: 

The  Lawrence  Livermore  National  Laboratory  Expert  Committee  has  reviewed  the 
methods  and  findings  of  the  site  investigations  conducted  at  George  Air  Force  Base 
(GAFB),  located  in  Victorville,  California,  with  regard  to  subsurface  releases  of 
aviation  fuel.  These  site  investigations  have  been  conducted  by  IT  Corporation  and 
others.  Included  is  our  evaluation  of  the  adequacy  of  the  site  conceptual  model  and 
risk  characterization  as  a  credible  basis  for  applying  risk-based  corrective  action 
(RBCA),  emphasizing  natural  attenuation  as  a  means  of  remediation,  to  achieve  site 
closure. 

This  letter  represents  the  first  of  two  deliverable  documents,  as  part  of  our  overall 
assessment.  It  is  intended  solely  as  a  brief  review  of  the  existing  site  data  and  risk 
characterization  models;  recommendations  regarding  additional  data  needed  to 
complete  the  site  model  review  are  provided.  Specifically,  the  review  included: 

•  An  assessment  of  site  characterization  and  the  contaminant  distribution 
model,  based  upon  conformance  with  standard  environmental  investigation 
practices  and  sound  scientific  judgment; 

•  A  preliminary  assessment  of  bioattenuation  signatures,  based  upon 
presentation  of  multiple  lines  of  evidence  such  as  natural  geochemical 
indicators  and  changes  in  plume  composition  over  time; 

•  An  assessment  of  risk  characterization,  based  upon  our  judgment  as  to  how 
significant  uncertainties  in  risks  posed  to  potential  receptors  have  been 
characterized. 

Detailed  analyses  of  key  assumptions  and  an  assessment  of  site-specific  natural 
attenuation  will  be  provided  in  a  follow-on  report  as  a  separate  deliverable  item. 

With  regard  to  GAFB,  the  dissolution  rate  of  the  free  product  lens  and  the  capacity 
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for  natural  attenuation  processes  to  prevent  further  downgradient  migration  of  the 
dissolved  hydrocarbon  plume  will  be  the  subject  of  focused  study. 

Site  Conceptual  Model 

In  our  opinion,  the  distribution  of  contaminants  at  GAFB  has  been  reasonably  well 
characterized.  Petroleum  hydrocarbons  (originally  JP-4  aviation  fuel)  are  present  in 
the  subsurface  environment  at  GAFB  as  a  result  of  leaks  from  underground  tanks 
and  piping  systems.  A  significant  volume  of  separate  phase  product,  estimated  to  be 
on  the  order  of  350,000  gallons,  persists  on  top  of  the  water  table,  covering  an  area  of 
approximately  32  acres,  with  product  thicknesses  ranging  from  a  thin  sheen  to 
several  feet.  A  plume  of  dissolved  fuel  hydrocarbon  compounds  (FHCs)  exists  in  the 
underlying  groundwater,  extending  a  few  hundred  feet  downgradient  of  the  free 
product  lens.  Isolated  pockets  of  low  concentrations  of  dissolved  halogenated 
hydrocarbons  (TCE,  PCE,  1,1-DCA)  have  also  been  noted  in  site  groundwater. 

The  site  hydrogeologic  model  is  well  developed  (37  monitoring  wells  have  been 
installed  onsite).  The  site  is  underlain  by  alluvial  sediments  consisting  of 
unconsolidated  sands,  silts,  clays,  and  gravels  in  an  arid  climatic  setting.  An 
unconfined  water-bearing  zone  underlies  the  site  (depth  to  the  water  surface  is  on 
the  order  of  100  feet);  this  water-bearing  zone  is  in  turn  underlain  by  a  clay  aquiclude 
and  a  regional  aquifer.  Hydraulic  conductivity  values  have  been  calculated  from  a 
number  of  pumping  tests  and  aquifer  tests  (although  analyses  associated  with  the 
two  methods  have  produced  estimated  values  which  are  not  entirely  self- 
consistent).  Taken  with  interpolated  hydraulic  head  values  and  porosity  estimates, 
these  data  suggest  groundwater  movement  in  a  north-northeasterly  direction  at  a 
rate  of  tens  to  hundreds  of  feet  per  year.  Little  recharge /evapotranspiration  is 
expected  in  association  with  the  unconfined  aquifer. 

The  spatial  extent  of  the  dissolved  FHC  constituents /BTEX  in  groundwater  appears 
to  reasonably  well-characterized  over  individual  sampling  rounds,  although  long¬ 
term  historical  trends  in  individual  wells  are  not  available.  Concentrations  of 
natural  attenuation  geochemical  indicator  species  (e.g.  dissolved  oxygen,  nitrate, 
sulfate,  iron,  methane,  redox  potential)  have  also  been  collected  from  16  monitoring 
wells  as  part  of  a  natural  attenuation  monitoring  study.  Clear  evidence  has  been 
presented  indicating  natural  bioattenuation  processes  are  occurring  at  GAFB.  These 
include  apparent  shrinkage  of  BTEX  constituent  plumes  with  time,  loss  of  BTEX 
components  in  comparison  to  recalcitrant  dissolved  FHC  constituents 
(trimethylbenzene  isomers),  and  changes  in  geochemical  indicator  species  which 
vary  in  space  in  a  manner  consistent  with  expected  bioattenuation  redox  regimes. 

Site  Risk  Characterization  and  Management  Strategies 

Standard  risk  assessment  analyses  were  conducted  for  the  GAFB  site.  Because  of  the 
depth  at  which  the  contamination  occurs,  the  present  land  use  at  GAFB,  and  the 
presumed  large  distances  to  potential  human  and  ecological  risk  receptors 


downgradient,  risks  to  human  health  through  common  exposure  pathways  and  to 
the  environment  appear  to  be  low.  However,  we  have  not  been  able  to  identify 
locations  and  screened  intervals  of  existing  water  supply  wells  downgradient  of  the 
GAFB  site  in  the  site  documentation.  Such  information  needs  to  be  provided  if 
available  to  apply  risk  based  corrective  action  to  the  site. 

Current  remediation  techniques  applied  at  GAFB  include  free  product  removal, 
capillary  fringe  bioventing,  and  natural  attenuation  of  dissolved  FHCs  in 
groundwater.  Free  product  remediation  consists  of  mobile  and  permanent  free 
product  recovery  systems,  which  remove  approximately  160  gallons  of  free  product 
per  day,  with  a  cumulative  total  of  approximately  32,000  gallons  to  date.  The 
bioventing  program  has  to  date  consisted  largely  of  a  treatability  test.  Natural 
attenuation  is  expected  to  provide  the  means  for  controlling  and  eventually 
remediating  the  dissolved  FHC  plume.  In  addition  to  presenting  parallel  lines  of 
evidence  supporting  the  occurrence  of  bioattenuation  at  GAFB,  modeling  studies 
using  BIOPLUME  II  were  conducted  by  IT  to  quantify  the  expected  decline  of  the 
hydrocarbon  plume  over  time  as  a  result  of  natural  attenuation.  Sensitivity  studies 
were  conducted  to  evaluate  the  significance  of  parameters  such  as  dispersivity  and 
anaerobic  decay  rate  on  estimated  time  to  achieve  cleanup  (i.e.  time  required  to 
reach  MCLs).  The  conclusion  of  these  modeling  exercises  was  that  cleanup  of  the 
dissolved  FHC  plume  would  be  achievable  within  40  years. 

In  principle,  the  application  the  BIOPLUME  II  model  for  evaluating  the  future 
behavior  of  the  dissolved  FHC  plume  at  GAFB  is  reasonable,  provided  that  (1) 
uncertainties  in  model  parameters  and  contaminant  sources  are  addressed,  and  (2) 
the  model  is  not  relied  upon  as  the  exclusive  means  for  assessing  future 
contaminant  behavior.  The  source  term  in  the  GAFB  BIOPLUME  II  model  is  based 
upon  quantifying  Fickian  diffusive  solute  fluxes  across  the  free  product-water 
interface.  These  calculations  are  very  sensitive  to  assumptions  regarding  the  mass 
of  the  free  product,  boundary  concentrations,  concentration  gradient,  and  so  forth. 
Our  preliminary  evaluation  of  the  free  product  dissolution  problem,  using  a  model 
based  on  a  dissolution  mixing  zone  factoring  in  regional  groundwater  flow,  suggests 
that  the  predicted  cleanup  time  projected  by  the  BIOPLUME  model  is  subject  to 
considerable  uncertainty.  In  our  opinion,  this  uncertainty  must  be  addressed  by 
additional  analyses,  as  the  estimated  time  that  natural  attenuation  can  be  expected  to 
remediate  the  site  depends  on  the  quantity  of  free  product  present  and  the  rate  of 
mass  removal.  In  response,  more  detailed  analyses  of  free  product  dissolution, 
using  alternative  models  (e.g.  time-dependent  diffusion,  mixing  zone)  will  be 
undertaken  by  the  Expert  Committee  for  inclusion  in  the  follow-on  report. 

Clarification  of  several  issues  concerning  the  free  product  lens  will  aid  in  our 
analyses.  Experience  suggests  that  free  product  removal  efforts  generally  do  not 
succeed  in  removing  more  than  about  30%  of  the  spill  volume.  Based  on  current 
free  product  estimates,  the  current  mass  removed  represents  approximately  10%  of 
the  estimated  free  product  source  mass.  Historical  free  product  removal  rates  may 
help  bound  the  uncertainty  in  the  source  mass.  Also,  we  note  that  the  estimates  of 


free  product  mass  present  in  the  subsurface  at  GAFB  have  tended  to  increase 
through  time,  ranging  from  roughly  250,000  gallons  to  well  over  400,000  gallons. 

We  recognize  that  this  may  be  a  result  of  further  characterization  as  the 
investigation  has  progressed.  However,  it  would  be  helpful  if  we  were  provided 
with  the  explanation  for  the  range  of  values  reported  would  also  be  provided  so  that 
we  may  better  understand  the  variability  of  the  free  product  mass.  Finally, 
compositional  variations  in  the  weathered  free  product,  in  comparison  to 
unweathered  JP-4,  may  provide  insights  into  the  rate  of  volatilization  of  lighter- 
weight  (and  therefore  more  soluble)  components. 

In  addition  to  addressing  the  source  term  uncertainty,  the  potential  for  natural 
attenuation  mechanisms  to  limit  the  extent  of  the  dissolved  FHC  plume  will  also  be 
evaluated  by  the  Expert  Committee  as  part  of  the  follow-on  site  report.  This 
evaluation  will  consist  of  applying  alternative  plume  bioattenuation  models  to  the 
GAFB  groundwater  data  set,  taking  into  account  geochemical  bioattenuation 
indicators.  The  results  of  these  analyses  will  be  used  in  an  attempt  to  validate  the 
natural  attenuation  conceptual  model. 

Summary 

Subsurface  investigations  conducted  at  GAFB  by  IT  Corporation  appear  to  be 
satisfactory  in  terms  of  characterizing  the  nature  of  the  subsurface  environment  and 
the  distribution  of  contaminants.  Thus,  we  do  not  foresee  any  need  to  collect 
additional  field  data  for  the  application  of  the  RBCA  decision-making  process.  The 
natural  attenuation  study  has  provided  credible  evidence  that  dissolved  FHCs  in 
groundwater  underlying  the  site  are  indeed  undergoing  bioattenuation  to  some 
extent.  Furthermore,  it  appears  that  existing  site  data  may  be  used  to  construct 
arguments  that  the  dissolved  FHC  plume  will  at  least  remain  stable  for  the 
foreseeable  future  (i.e.  the  plume  in  all  likelihood  will  not  advance  further  in  the 
downgradient  direction).  This  scenario  will  be  investigated  by  the  Expert 
Committee  using  additional  plume  analyses,  including  a  quantitative  assessment  of 
bioattenuation  signatures. 

In  our  opinion,  the  cleanup  time  estimates  provided  by  modeling  analyses  using 
BIOPLUME  II  are  characterized  by  significant  uncertainty.  This  is  primarily  because 
of  potential  uncertainties  with  regard  to  the  source  term  associated  with  the  large 
free  product  volume.  As  a  result,  given  the  existing  data,  it  is  difficult  to  ascertain 
the  length  of  time  shallow  groundwater  immediately  underlying  the  site  will 
continue  to  be  impacted  by  petroleum  hydrocarbons.  The  Expert  Committee  will 
address  the  uncertainty  associated  with  the  time  required  for  free  product 
dissolution  by  additional  quantitative  analyses  in  the  follow-on  site  report.  To  assist 
in  this  effort,  the  following  data  are  requested: 

(1)  Historical  free  product  recovery  estimates; 

(2)  A  brief  summary  explanation  of  the  uncertainty  in  free  product  mass 
estimates; 


(3)  Chromatographs  from  any  historical  free  product  analysis  to  evaluate  product 
weathering. 

Also,  to  facilitate  risk  analyses  as  part  of  the  RBCA  decision-making  process,  we 
request: 

(4)  Downgradient  water  supply  well  locations  and  screened  intervals  (if 
available). 

We  look  forward  to  preparing  the  GAFB  site  report  following  receipt  of  this 
information. 


Sincerely, 


Walt  W.  McNab,  Jr. 

Expert  Committee  Lead  for  George  A.F.B. 


cc:  John  Adams,  SWRCB 
Kandi  Brown,  IT  Corp 
Jehiel  Cass,  Lahontan  RWQCB 
James  Chang,  US  EPA  Region  IX 
Steve  Cullen,  UCSB 
Lome  Everett,  UCSB 
Marty  Faile,  USAFCEE 
James  Giannopoulos,  SWRCB 
Juan  Jimenez,  DTSC 
Paul  Johnson,  ASU 
William  Kastenberg,  UCB 
Michael  Kavanaugh,  Environ 
Charles  Lee,  USAFCEE 
Richard  Ragaini,  LLNL 
Matt  Small,  US  EPA  Region  IX 
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Free  Product  Dissolution/Volatilization  Analysis 

B-l.  Background 

The  source  term  used  to  describe  the  contribution  of  BTEX  components  from  dissolution  of  the 
non-aqueous  phase  liquid  (NAPL)  pool  in  the  BIOPLUME  II  model  involved  the  calculation  of 
steady-state  diffusive  fluxes  based  on  Fick’s  law  and  then  emulating  the  fluxes  with  a  series  of 
injection  wells  operating  at  very  low  flow  rates.  Using  an  assumed  concentration  gradient  across 
the  NAPL-groundwater  interface,  diffusive  fluxes  were  assumed  to  persist  at  constant  levels  until 
BTEX  components  were  depleted.  The  assumption  of  steady-state  influxes  of  BTEX  components 
into  the  dissolved  phase  is  only  a  crude  approximation  because  concentration  gradients,  and  hence 
solute  fluxes,  across  the  interface  will  tend  to  decrease  over  time,  but  persisting  for  a  longer  period 
than  predicted  by  a  steady-state  model.  Nevertheless,  the  assumption  of  a  steady-state  flux 
provides  a  reasonable  order-of-magnitude  approach  at  quantifying  the  problem.  However,  the 
specific  approach  used  in  the  BIOPLUME  II  model  to  simulate  diffusive  fluxes  of  BTEX  across 
the  NAPL-groundwater  boundary  is  problematic  for  two  principle  reasons: 

1 .  The  definition  of  concentration  gradient  is  arbitrary.  Fluxes  are  calculated  based  upon  a 
concentration  gradient  defined  over  a  distance  of  20  ft  between  the  concentrations  of  BTEX 
components  in  the  JP-4  NAPL  phase  and  the  highest  concentrations  of  BTEX  measured  in 
groundwater  monitoring  wells. 

2.  Limited  solubility  of  BTEX  components  in  the  aqueous  phase  is  not  considered.  This  is  a 
key  omission;  pure-phase  BTEX  components  are  characterized  by  solubilities  in  water  that 
are  less  than  concentrations  in  JP-4  NAPL  (e.g.,  3,880  mg/L  benzene  concentration  in 
NAPL — based  on  a  weight  fraction  of  0.5  percent — versus  a  pure-phase  solubility  of 
1,770  mg/L).  In  addition,  because  JP-4  NAPL  is  a  mixture  of  a  large  number  of 
components,  solubilities  of  individual  species  such  as  benzene  will  be  less  than  those  of 
their  pure-phase  equivalents  as  approximated  by  Raoulf  s  law1,  ignoring  cosolvency  effects 
(see  discussion  below).  As  a  result,  concentration  gradients  calculated  using  BTEX 
concentrations  in  the  NAPL  phase  are  not  physically  meaningful. 


‘Raoult’s  law  states  that  the  partial  pressure  of  a  given  species  in  the  gas  phase  is  equal  to  the  mole  fraction  of  the 
component  in  the  liquid  phase  multiplied  by  the  partial  pressure  of  the  pure  phase  substance.  By  analogy,  the 
aqueous  solubility  of  a  component  in  a  NAPL  mixture  (e.g.,  benzene)  is  equal  to  its  mole  fraction  in  the  liquid 
phase  multiplied  by  the  aqueous  solubility  of  the  pure  substance.  This  idealization  does  ignore  cosolvency  effects, 
which  tend  to  raise  the  solubility  of  a  component  above  that  predicted  by  Raoult’s  law.  However,  in  the  case  of  JP- 
4,  it  is  likely  that  benzene  and  the  other  BTEX  components  are  not  subject  to  strong  cosolvency  effects,  given  that 
highly  soluble  additives  such  as  methanol  or  methyl  tertiary  butyl  ether  (MTBE)  are  not  present  in  the  mixture. 
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B-2.  Plume  Dissolution  Model 

A  more  rigorous  approach  is  to  recognize  that  1)  aqueous  phase  solubilities  of  individual 
components  in  the  NAPL  mixture  control  diffusive  transport  across  the  NAPL-water  interface,  and 
2)  concentration  gradients  across  the  NAPL-water  interface  are  influenced  by  the  lateral  advective 
transport  of  dissolved  hydrocarbons  from  dissolution  of  any  NAPL  present  in  the  upgradient 
direction.  Hunt  et  al.  (1988)  presented  a  steady-state  conceptual  model  of  the  dissolution  of  a 
NAPL  pool  as  described  by  the  equation, 

dc  d2c 

V  ~rh  ~  ^72 

dx  dz  (B-l) 

subject  to  the  boundary  conditions, 

c(0  <  x  <  L,  z  =  0)  =  cs 

c(x,  z  =  °o)  =  0 

c(x  =  0,  z)  =  0 

where  v  refers  to  the  groundwater  pore  velocity,  dc/dx  the  concentration  gradient  in  the  horizontal 
direction,  Dz  the  effective  dispersion  coefficient  in  the  vertical  direction,  and  d2c/dz2  the  gradient  of 
the  concentration  gradient  in  the  vertical  direction.  Equation  B-l  states  that  the  downward 
dispersive  vertical  movement  of  the  dissolved  contaminant  at  any  point  below  the  NAPL-water 
interface  is  equal  to  the  lateral  advective  flux.  The  solution  to  this  equation  is: 


c(L,z)  =  Cserfc 


(  DZL\ 


L  V  v  ; 


1/2 


(B-2) 


where  L  refers  to  the  distance  along  the  NAPL  product  pool  from  the  upgradient  end,  z  the  vertical 
depth  below  the  NAPL-water  interface,  and  Cs  the  effective  aqueous  solubility  of  the  NAPL 
component.  Equation  B-2  may  be  used  to  evaluate  solute  fluxes  across  the  NAPL-water  interface 
by  numerically  differentiating  concentration  with  distance  to  obtain  a  concentration  gradient 
estimate  at  point  L  along  the  base  of  the  NAPL  plume: 


dc  _  c(L,&)-c(L,0) 
dz j=o  Sz 


(B-3) 


where  8z  is  some  small  distance  (e.g.,  1  x  10  5  m),  and  then  multiplying  by  the  effective 
dispersion  coefficient: 


mass  _  dc 
time  z  dzz= o 


(B-4) 
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where  A  refers  to  a  horizontal  cross-sectional  area.  Dissolution  time  can  then  be  calculated  by 
assuming  a  mean  thickness  for  the  NAPL  pool  and  hence  the  benzene  mass  per  unit  horizontal 
cross-sectional  area. 

Table  B-l.  Parameters  for  dissolution  model  of  Hunt  et  al.  (1988). 


Parameter 


Value 


Notes 


Benzene  mole  fraction  in  JP-4 


Dz  (effective  vertical  dispersion 
coefficient) 

K,  hydraulic  conductivity 

Vh,  hydraulic  gradient 
<(>,  porosity 

A  (area  of  NAPL  pool) 

V  (volume  of  NAPL  pool) 


Lognormal.  Mean  =  0.56%; 
standard  deviation  =  0.28%. 


Lognormal.  Mean  =  4.4  x  10"9 
m2/sec;  standard  deviation  = 

8.91  x  10“9  m2/sec.  Corresponds 
to  the  lst-percentile  =  10"10 
m2/sec  and  the  90th  percentile  = 
lCf8  m2/sec. 

Lognormal.  Mean  =  61.2  ft/day; 
standard  deviation  =  30.4 
ft/day.  Corresponds  to  the  10th- 
percentile  =  30  ft/day  and  the 
90th-percentile  =  100  ft/day. 

Lognormal.  25th-percentile  = 
0.002;  75th-percentile  =  0.003. 

Normal.  Mean  =  0.25;  standard 
deviation  =  0.03. 

Normal  distribution.  Mean  = 
1,860,000  ft2;  standard  deviation 
=  186,000  ft2  (assumed  to  be  equal 
to  10%  of  the  mean  for 
simplicity). 

Normal  distribution.  Mean 
=350,000  gal;  standard  deviation 
=  35,000  gal  (assumed  to  be  equal 
to  10%  of  the  mean  for 
simplicity). 


Pure  benzene  solubility  =  1770 
mg/L;  Raoult's  law 
approximation;  mole  fraction  of 
benzene  assumed  equal  to  its 
weight  fraction  in  JP-4,  0.5%  (IT, 
1996).  Variations  in  JP-4 
composition  and  departures  from 
the  ideality  of  Raoult's  law 
may  explain  slightly  higher 
concentrations  of  benzene  noted 
in  some  GAFB  groundwater 
samples. 

Based  on  vertical  mixing  model 
and  Peclet  number 
considerations.  See  discussion  in 
text. 


Based  upon  site  data  (IT,  1996). 


Based  upon  site  data  (IT,  1996). 
Postulated. 

IT  (1996). 


IT  (1996). 


A  Monte  Carlo  simulation  was  conducted  using  Equations  B-2  through  B-4  to  provide  insight 
into  the  amount  of  time  required  for  benzene  to  dissolve  out  of  the  free  product  pool  at  GAFB  by 
assuming  the  system  remains  at  steady-state.  Probability  distributions  for  the  governing  parameter 
are  shown  on  Table  B-l;  values  were  estimated  from  site  documentation  (e.g.,  IT,  1992)  or 
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otherwise  estimated  based  on  judgment  of  likely  site  conditions.  Modeling  consisted  of 
programming  Equations  B-2  through  B-4  into  a  Microsoft  Excel  spreadsheet  and  then  using  the 
Crystal  Ball  add-in  package  (Decisioneering,  Inc.,  1996)  for  Excel  to  run  a  total  of  10,000 
realizations  with  the  prescribed  parameter  distributions. 

Effective  solubility  of  benzene,  Cs,  was  calculated  based  on  the  pure  phase  solubility  of 
benzene,  Cp,  (1,770  mg/L)  and  a  Raoult’s  law  approximation  using  the  mole  fraction  of  benzene  in 
JP-4,  Xb, 

Cs  =  CpXb  (B-5) 

Groundwater  velocity  was  calculated  using  Darcy’s  law, 

fCVh 


where  K  is  the  hydraulic  conductivity,  V/t  the  hydraulic  conductivity,  and  0  the  porosity. 
Assuming  a  circular  NAPL  pool  as  a  simplifying  assumption,  a  representative  distance  along  the 
base  of  the  pool  representing  an  average  dissolution  rate  can  be  assumed  to  equal  the  radius  of  the 
pool,  thus  the  parameter  L  in  Equation  (B-2)  was  estimated  by: 


L  = 


(B-7) 


where  A  is  the  area  of  the  NAPL  pool  in  plan  view. 

Cleanup  times  were  forecast  by  dividing  the  benzene  mass  per  unit  area  in  the  NAPL  pool  by 
the  modeled  flux  rate  (as  calculated  by  Equation  B-4).  Benzene  mass  per  unit  area  was  estimated 
by  first  dividing  the  NAPL  volume  by  its  area  to  yield  the  NAPL  pool  mean  thickness,  then 
assuming  a  density  of  0.8  g/cm3  to  convert  to  mass  per  unit  area,  and  finally  correcting  this  value 
by  the  benzene  mole  fraction  (assumed  equivalent  to  the  mass  fraction).  The  cumulative 
distribution  function  resulting  from  the  Monte  Carlo  analyses  for  total  dissolution  time  of  benzene 
out  of  the  NAPL  phase  is  shown  on  Figure  B-l.  The  median  value  is  on  the  order  of  50  years. 
Because  of  the  steady-state  assumption,  the  estimates  produced  are  conservative  because  benzene 
concentrations  in  the  NAPL,  and  hence  along  the  saturation  boundary  layer,  in  reality  would 
decline  with  time.  This  would  result  in  a  decrease  in  gradient  and  hence  a  reduction  in  the 
dissolution  rate. 


B-2.1.  Effective  Vertical  Dispersion  Coefficient 

Vertical  fluxes  calculated  by  the  methodology  outlined  above  show  a  high  degree  of  sensitivity2 
to  the  vertical  dispersion  coefficient,  Dz,  shown  in  Equation  B-l.  In  principle,  this  value  could 


Sensitivity  analysis  of  the  Monte  Carlo  simulation  output,  based  on  contribution  to  variance  in  rank  correlation, 
suggests  that  the  effective  vertical  dispersion  coefficient  has  the  most  significant  effect  by  far  on  cleanup  time 
estimates  (approximately  85%  of  the  variance  is  attributable  to  uncertainty  in  Dz),  followed  by  hydraulic  conductivity 
(11%),  NAPL  volume  ( 1 .5%),  hydraulic  gradient  (1  %),  and  NAPL  area  ( 1  %). 
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range  from  the  effective  molecular  diffusion  coefficient,  on  the  order  of  1  x  1CT10  m2/sec,  to  several 
orders-of-magnitude  higher  if  vertical  hydrodynamic  dispersion  is  a  significant  process  at  the 
NAPL-water  interface.  One  method  for  evaluating  a  realistic  range  of  possible  vertical  dispersion 
coefficients  is  to  examine  the  predicted  concentrations  in  the  underlying  aqueous  phase  as  a 
function  of  depth  according  to  Equation  B-2.  Because  groundwater  samples  are  not  usually 
collected  at  discrete  points  but  rather  are  a  reflection  of  integrated  measurements  over  the  length  of 
a  well  screen,  it  is  appropriate  to  consider  depth-averaged  concentrations  which  could  be  measured 
under  the  steady-state  assumptions  used  in  the  development  of  Equation  B-2. 

Calculated  vertically-averaged  concentrations  of  benzene  (at  plume  center)  along  discrete  depth 
intervals  as  a  function  of  the  vertical  dispersion  coefficient  are  shown  on  Table  B-2  (assuming  Cs  = 
8.85  mg/L,  K  =  60  ft/day,  V/i  =  0.002,  <)>  =  0.25,  and  L  =  230  m).  Given  historical  measurements 
of  benzene  concentrations  in  groundwater  at  GAFB  provided  in  site  documentation,  it  appears  that 
values  of  Dz  greater  than  1  x  10  7  m2/sec  produce  vertically-averaged  concentrations  of  benzene 
that  exceed  observation  to  a  considerable  degree,  especially  when  averaged  over  a  depth  of  10  m. 


Table  B-2.  Vertically-averaged  concentrations  parts  per  billion  (ppb)  versus  effective 
vertical  dispersion  coefficient. 


Depth 

(m) 

1  x  10'10 
m2/sec 

1  x  10'9 
m2/sec 

1  x  10-8 
m2/sec 

1  x  10"7 
m2/sec 

0 

8,850 

8,850 

8,850 

8,850 

0.1 

6,850 

8,183 

8,638 

8,783 

0.5 

2,755 

5,769 

7,799 

8,515 

1 

1,503 

3,719 

6,796 

8,181 

2 

787 

1,986 

5,072 

7,525 

5 

324 

818 

2,390 

5,731 

10 

164 

413 

1,208 

3,028 

A  second  approach  at  quantifying  Dz  is  to  consider  the  nature  of  the  mixing  regime.  This  is 
accomplished  through  evaluating  the  dimensionless  Peclet  number,  which  is  defined  as  average 
pore  velocity  multiplied  by  average  throat  velocity  divided  by  the  effective  molecular  diffusion 
coefficient.  The  effective  molecular  diffusion  coefficient  is  defined  as  molecular  diffusion 
coefficient  divided  by  a  tortuosity  factor.  The  Peclet  number  is  a  measure  of  the  contribution  to 
mixing  of  mechanical  mixing  compared  to  molecular  diffusion.  At  values  much  less  than  1.0, 
molecular  diffusion  acts  as  the  dominant  mixing  mechanism  and  mechanical  dispersion  may  be 
neglected.  At  values  much  greater  than  1.0  the  opposite  is  true.  When  the  Peclet  number  is  in  the 
range  of  0.5  to  5  it  must  be  assumed  that  both  molecular  diffusion  and  mechanical  dispersion  are 
significant  and  must  be  accounted  for.  Various  Peclet  numbers  for  a  range  of  throat  diameters  are 
listed  on  Table  B-3  and  B-4  using  site  specific  parameters  for  the  GAFB  site  and  the  properties  of 
benzene. 
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Table  B-3.  Peclet  number  parameters. 


Molecular  diffusion  coefficient 
(m2/s) 

8.3  x  10"10 
m2/sec 

Tortuosity 

2 

Porosity 

0.2 

Effective  diffusion  coefficient  (cm2) 

8.3  x  10~n  m2/sec 

Average  velocity  (cm/s) 

2.5  x  10  ~7  m/sec 

Table  B-4.  Peclet  number  versus  pore  diameter. 


Pore  diameter  (m) 

Peclet  No. 

1  x  10'3 

3.0 

1  x  10-4 

0.3 

1  x  10-5 

0.03 

These  results  indicate  that  mixing  is  probably  dominated  by  molecular  diffusion.  In  the 
unlikely  case  that  the  throat  diameter  is  as  large  as  1 .0  millimeter,  the  Peclet  number  ranges  into  the 
region  where  mechanical  mixing  should  be  considered  along  with  diffusion.  However,  at  more 
realistic  throat  diameters  of  1/10  to  1/100  of  a  millimeter,  mixing  is  completely  dominated  by 
molecular  diffusion  and  mechanical  mixing  can  be  neglected.  Thus,  taken  with  dissolved  BTEX 
concentrations  measured  in  groundwater  near  the  NAPL  plume,  these  analyses,  suggest  that  the 
vertical  dispersion  coefficient  is  generally  small.  For  modeling  purposes,  therefore,  a  lognormal 
distribution  of  the  effective  vertical  dispersion  coefficient  between  1  x  10-8  and  1  x  10-1^m2/sec 
was  assumed. 

B-2.2.  An  Alternative  Dissolution  Model  Based  Upon  a  Mixing  Zone 
Assumption 

As  an  alternative  means  for  estimating  NAPL  dissolution  time,  a  second  modeling  approach, 
based  on  a  mixing  zone  concept  for  the  NAPL-groundwater  interface,  was  also  evaluated.  This 
approach  simply  assumes  that  a  finite  thickness  mixing  zone  exists  in  groundwater  immediately 
below  the  NAPL  pool,  which  is  saturated  with  benzene3.  In  this  case,  dissolved  contaminants  are 
flushed  away  from  the  NAPL  pool  at  a  rate  equal  to  the  concentration,  times  the  advective  pore 
velocity,  times  the  cross-sectional  area  of  the  mixing  zone  underneath  the  plume.  A  Monte  Carlo 
analyses  was  performed  using  the  same  constraints  as  the  first  model  for  gradient,  porosity, 
hydraulic  conductivity,  and  NAPL  pool  volume  and  area.  Two  probability  distribution  models 
were  used  to  described  the  thickness  of  the  mixing  zone.  The  first  model  (Case  1)  assumed  an 


Both  models  are  insensitive  to  the  weight  fraction/mole  fraction  of  benzene,  assuming  a  Raoult's  law  relationship. 
For  example,  if  the  weight  fraction  of  benzene  in  JP-4  were  2%  instead  of  0.5%,  its  solubility,  according  to  Raoult’s 
law,  would  be  35.4  mg/L,  so  the  flux  out  the  NAPL  phase  would  be  four  times  as  high.  At  the  same  time, 
however,  there  would  be  four  times  as  much  benzene  to  remove,  so  the  effects  of  higher  benzene  concentrations 
would  tend  to  cancel  one  another. 
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exponential  distribution  of  mixing  zone  thickness,  with  a  rate  coefficient  of  2.0.  The  second  model 
(Case  2)  assumed  a  uniform  probability  distribution  for  the  mixing  zone  thickness,  ranging  from 
0.01  m  to  1.0  m.  The  results  of  both  of  these  simple  models,  illustrated  on  Figure  B-l,  also 
suggest  a  long  dissolution  time,  generally  much  greater  than  100  years. 

As  expected,  sensitivity  analyses  reveal  that  the  uncertainty  regarding  the  depth  of  the  mixing 
zone  is  the  primary  factor  in  influencing  variance  in  cleanup  time  estimates  (84%  and  69%  of  the 
contribution  to  variance  for  Case  1  and  Case  2,  respectively)  followed  by  much  smaller 
contributions  due  to  uncertainties  in  hydraulic  conductivity,  gradient,  and  NAPL  pool  volume. 

B-2.3.  Volatile  Losses  Into  the  Vadose  Zone 

In  environments  not  characterized  by  large  recharge,  the  transport  of  chemicals  in  the  vadose 
zone  is  mainly  a  diffusive  process.  Mayer  et  al.  (1974),  employing  a  methodology  by  which  water 
transport  was  assumed  to  be  negligible,  proposed  five  different  models  for  pesticide  movement  due 
to  diffusive  forces  alone.  The  methodology  used  a  separation  of  variables  approach  as  discussed 
by  Carslaw  and  Jaeger  (1959)  and  has  been  used  by  others  as  a  screening  level  tool  (Jury  et  al., 
1983,  1984a,b,c,  and  1990).  A  similar  approach  using  boundary  conditions  appropriate  for  the 
site  is  discussed  below. 

Vapor  phase  transport  of  benzene  from  the  NAPL  plume  was  examined  to  determine  the  time 
frame  in  which  it  would  be  removed.  The  equation  of  diffusive  transport,  written  in  terms  of  the 
total  concentration  is: 


dQ  _  £>  <?2Cr 

*  dz2 

with  the  total  concentration  being: 


CT 


=  C„ 


0. 


M±] 
H  J 


(B-8) 


(B-9) 


01  is  the  volumetric  water  content4,^  is  the  vapor  content,  H  is  the  dimensionless  Henry’s  law 
constant,  pt,  is  the  bulk  soil  density,  and  Kj  is  the  distribution  coefficient  for  the  chemical  (based 
on  a  reported  value  by  IT  Corp.,  1996). 

The  effective  diffusion  coefficient  is  based  on  the  gas  diffusion  coefficient  (the  water  diffusion 
coefficient  has  been  neglected),  and  is: 


4Sediments  in  the  vadose  zone  environment  at  GAFB  consist  of  silt  and  sand,  with  an  average  effective  porosity  of 
0.18.  The  stated  volumetric  water  content  provided  in  the  GAFB  Operable  Unit  2  Natural  Attenuation  Program 
(LLNL  Expert  Committee),  May  1996  report  is  a  misnomer  (Larry  Tyner,  Personal  communication,  IT  Corp., 
1996).  The  volumetric  water  content  is  given  as  an  average  0.158,  distributed  from  0.15  to  0.23.  Based  on  the 
given  ranges  of  porosity,  this  would  correlate  to  a  very  moist  vadose  zone  environment,  in  an  area  characterized  by 
low  rainfall  rates  and  a  large  concrete  cap  covering  most  of  the  site.  It  is  therefore  assumed  that  specified  water 
content  is  actually  the  saturated  fraction,  resulting  in  a  range  of  volumetric  water  content  from  0.023  to  0.078. 
Questions  remain  as  to  the  methodology  of  the  tests  performed  to  determine  saturated  fraction  and  porosity. 
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Deff=  — 

eff  Rr 


(B-10) 


The  unsaturated  retardation  factor,  R,  is  defined  by: 

r  =  1  +  Jl.+B£l 

He,  He, 


(B- 11) 


and  the  tortuosity,  x,  for  gaseous  diffusion  in  air-filled  pore  spaces  only,  is  the  Millington  (1959) 
equation  corrected  by  Falta  et  al.  (1992), 


T=i_ 


(B-12) 


The  boundary  conditions  for  the  problem  is  one  of  assumed  zero  flux  across  the 
low-permeability  concrete  cap,  with  a  constant  concentration  at  the  NAPL-vadose  zone  interface. 
The  NAPL  dimensional  area  is  large  enough  in  comparison  to  the  vadose  zone  thickness  that  a 
one-dimensional  (1-D)  idealization  is  for  estimating  a  time  scale  necessary  to  reach  steady  state  in 
the  vadose  zone.  The  boundary  conditions  are  shown  in  the  conceptual  scenario  (Fig.  B-2), 


/  =  0  CTo=  0  0  <z<L 

x  =  0  CTi  =  q 

x  =  L  _  o 

oz 

Using  a  separation  of  variables  approach,  it  may  be  shown  that 
C-CTi  ^  2(-l )" 


-(n  +  0.5)V 


LVA 


rff 


COS 


(«  +  0.5);r- 


(B-13) 


The  result  of  this  series  is  shown  in  Figure  B-3.  This  equation  has  been  dimensionalized,  so 
that  the  figure  is  representative  of  any  concentration  at  x  =  0,  or  for  any  uniform  concentration  that 
may  be  present  within  the  vadose  zone.  For  this  problem,  the  initial  concentration  within  the 
vadose  zone  is  taken  to  be  zero.  For  a  volumetric  water  content  of  0.028,  after  approximately 
70  years,  the  concentration  is  within  95%  of  the  steady-state  condition  possible  within  the  vadose 
zone  above  the  plume.  The  equivalent  steady  state  is  315  years  for  a  water  content  of  0.078. 
These  steady  state  conditions  give  a  gas  phase  concentration  for  benzene  of  approximately 
2.5  parts  per  million  (ppm). 

The  vapor  phase  concentration  of  2.5  ppm,  which  is  the  driving  force,  is  derived  from  an 
assumption  of  a  constant  surface  concentration  at  the  NAPL  surface.  Although  this  is  recognized 
to  be  conservative,  this  estimation  is  primarily  performed  to  develop  an  “order  of  magnitude” 
model  of  the  relative  contributions  of  dissolution  vs.  volatilization. 
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The  total  contaminant  mass  that  the  vadose  zone  has  the  capacity  to  assimilate  is  derived  from 
the  total  concentration  times  the  volume.  It  is  found  in  the  above  conservative  estimate  that  the 
vadose  zone,  idealized  as  a  one-dimensional  slab,  has  the  capacity  to  absorb  approximately 
1.4  times  the  benzene  assumed  to  be  resident  within  the  NAPL.  Using  this  constant  driver 
concentration,  the  vadose  zone  is  capable  of  absorbing  the  full  amount  of  benzene  in  25  years  for 
the  water  content  of  0.028,  and  100  years  for  the  water  content  of  0.078. 

This  analysis  does  not  take  account  into  several  factors,  including  constituents  already  within 
the  vadose  zone,  radial  transport  outward  at  the  fringes,  variable  water  content,  or  possible 
degradation.  Due  to  the  relatively  impermeable  boundary  condition  at  the  surface,  and  therefore  a 
low  diffusive  flux  of  oxygen,  degradation  may  be  oxygen  limited  near  the  center  of  the  plume, 
varying  outwards  towards  the  edges.  There  are  10  biovents  that  are  installed  around  the  site  above 
the  plume  that  are  providing  oxygen,  but  without  design  details  of  these  biovents,  the  oxygen 
content  within  the  vadose  zone  is  currently  unknown. 
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Figure  B-1.  Forecast  cumulative  distributions  of  free  product  dissolution  time. 
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Conceptual  scenario  of  volatilization/dissolution  of  NAPL  plume. 
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Appendix  C 


Plume  Natural  Attenuation  Analysis 

C-l.  Overview 

As  part  of  a  natural  attenuation  study  conducted  at  GAFB,  a  number  of  potential  indicators  of 
bioattenuation  of  the  dissolved  hydrocarbons  have  been  evaluated  during  site  investigation 
activities.  Measured  variations  in  dissolved  concentrations  of  hydrocarbon  constituents  are  not 
inconsistent  with  significant  bioattenuation.  However,  it  is  very  difficult  to  conclusively  assert  that 
the  plume  movement  is  being  controlled  by  bioattenuation  based  solely  on  measured  changes  in 
concentration,  because  the  sampling  history  does  not  provide  sufficient  temporal  and  spatial 
resolution  over  a  sufficiently  long  period  of  time.  In  this  case,  indirect  methods  of  evaluating 
plume  stability,  using  geochemical  indicators  and  relatively  simple  analytical  solute  transport 
models,  may  be  employed  to  provide  insight  into  plume  stability. 

C-2.  Plume  Length  Analysis 

Rice  et  al.  (1995)  noted  that  dissolved  benzene  typically  stabilized  in  most  groundwater 
environments,  with  the  length  of  benzene  plumes  ranging  from  100  ft.  (50th  percentile)  to  260  ft. 
(90th  percentile)  among  the  sites  examined.  This  plume  length  stabilization  reflects  a  balance 
between  dissolved  contaminant  mass  introduced  from  residual  sources  and  attenuation  processes 
(particularly  biodegradation)  integrated  across  the  plume. 

To  address  the  issue  of  plume  stability  in  the  special  case  of  GAFB,  Monte  Carlo  analyses 
were  performed  on  forecast  plume  lengths  using  the  analytical  solution  to  the  advective-dispersive 
solute  transport  equation  provided  by  Domenico  (1987).  The  two-dimensional  form  of  this 
solution  is, 
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where  C0  refers  to  the  source  concentration,  ax  the  longitudinal  dispersivity,  ay  the  transverse 
dispersivity ,  A.  the  first-order  decay  coefficient,  v  the  groundwater  pore  velocity,  R  the  retardation 
coefficient,  Y  the  width  of  the  line  source,  t  the  monitor  time,  and  x  and  y  the  coordinates  of  the 
monitor  point  relative  to  the  source.  The  model  used  in  this  evaluation  consisted  of  programming 
Equation  (C-l)  into  a  Microsoft  Excel  spreadsheet  and  then  using  the  Ciystal  Ball  add-in  package 
(Decisioneering,  Inc.,  1996)  for  Excel  to  run  Monte  Carlo  simulations  using  probability 
distributions  associated  with  the  given  parameters  (Table  C-l). 

Forecasting  probable  plume  lengths  associated  with  a  continuous  source  is  achieved  by 
predicting  concentrations  downgradient  of  the  source  (centerline)  and  then  using  a  suitable  search 
algorithm  such  as  bisection  or  Newton’s  method  to  find  the  location  downgradient  at  which  a 
certain  concentration  is  reached.  Plume  length  may  then  be  defined  as  the  distance  from  the  source 
to  such  a  location.  This  technique  was  used  to  generate  forecast  cumulative  distributions  of 
benzene  plume  lengths  as  defined  by  the  10  ppb  contour. 

The  constant  concentration  at  the  line  source  is  based  on  the  effective  solubility  of  benzene,  C5, 
which  was  calculated  using  the  pure  phase  solubility  of  benzene,  C  ,  (1,770  mg/L)  and  Raoult’s 
law  approximation  using  the  mole  fraction  of  benzene  in  JP-4,  Xb, 

c0  =CS  =  CpXb  (C-2) 

For  simplicity,  the  line  source  width  term  in  Equation  (C-l)  is  assumed  to  reflect  the  width  of  the 
observed  NAPL  plume  (idealized  as  a  square),  where  the  layer  of  water  in  immediate  contact  with 
the  NAPL  is  assumed  to  be  saturated  with  benzene  to  its  effective  solubility  limit.  Thus, 


Y  =  jA  (C-3) 

where  A  is  the  area  of  the  NAPL  pool.  The  retardation  coefficient,  R,  is  calculated  through  the 
expression, 


i?=l  + 


X  ocfocP 
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(C-4) 


where  Koc  is  organic  carbon  partitioning  coefficient  (83  mL/g  for  benzene),  foc  the  organic  carbon 
fraction  of  the  sediment,  pb  the  soil  bulk  density,  and  <j>  the  porosity. 

Flow  velocity  is  calculated  for  each  realization  in  accordance  with  Darcy’s  law, 

KXh 


In  studies  conducted  at  GAFB  with  an  in  situ  flow  meter  device,  Wilson  et  al.  (1995)  found 
that  local  preferential  flow  pathways  existed  within  the  hydrocarbon-impacted  zone  of  the  aquifer, 
characterized  by  flow  velocities  an  order-of-magnitude  higher  than  those  calculated  using  Darcy’s 
law  and  hydraulic  conductivity  estimates  based  upon  slug  tests  and  pumping  tests.  However, 
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Table  C-l.  Assumptions  used  in  probability  forecast  model. 


Parameter 

Benzene  mole  fraction  in  JP-4 


K,  hydraulic  conductivity 


Vh,  hydraulic  gradient 
<t>,  porosity 

foc,  fractional  organic  carbon 
content 


X,  decay  coefficient 


ax  ratio  (ratio  of  longitudinal 
dispersivity  to  plume  length 
scale,  vt ) 

ay  ratio  (ratio  of  transverse 
dispersivity  to  longitudinal 
dispersivity) 

t  (time  elapsed  between 
introduction  of  source  and 
groundwater  quality  sample) 


A  (area  of  NAPL  pool) 


Probability  distribution 

Lognormal.  Mean  =  0.56%; 
standard  deviation  =  0.28%. 


Lognormal.  Mean  =  61.2  ft/day; 
standard  deviation  =  30.4 
ft/day.  Corresponds  to  the  10th- 
percentile  =  30  ft/day  and  the 
90th-percentile  =  100  ft/day. 

Lognormal.  25th  percentile  = 
0.002;  75th-percentile  =  0.003. 

Normal.  Mean  =  0.25;  standard 
deviation  =  0.03. 

Lognormal.  Mean  =  0.08%; 
standard  deviation  =  0.29%. 
Corresponds  to  the  5th- 
percentile  =  0.02%  and  the  95th 
percentile  =  0.29%. 

Lognormal.  Mean  =  0.86%  day"1; 
standard  deviation  =  2.0%  day"1. 


Lognormal.  Mean  =  0.13; 
standard  deviation  =  0.10 

(median  =  0.10). 

Lognormal.  Mean  =  0.13; 
standard  deviation  =  0.10 

(median  =  0.10). 

Normal.  Mean  =  27  years  (10,000 
days);  standard  deviation  =  5.4 
years(20%  of  mean;  truncated  at 
30  years).  Doubled  for  projection 
in  year  2022. 

Normal  distribution.  Mean  = 
1,860,000  ft2;  standard  deviation 
=  186,000  ft2  (assumed  to  be  equal 
to  10%  of  the  mean  for 
simplicity). 


Basis 

Pine  benzene  solubility  =  1770 
mg/L;  Raoult's  law 
approximation;  mole  fraction  of 
benzene  assumed  equal  to  its 
weight  fraction  in  JP-4,  0.5%  (IT, 
1996).  Variations  in  JP-4 
composition  and  departures  from 
the  ideality  of  Raoult's  law 
may  explain  slightly  higher 
concentrations  of  benzene  noted 
in  some  GAFB  groundwater 
samples. 

Based  upon  site  data  (IT,  1996). 


Based  upon  site  data  (IT,  1996). 
Postulated. 

Based  upon  site  date  (IT,  1996). 


Based  on  decay  coefficients 
suggested  by  Wilson  et  al.  (1995) 
for  GAFB.  Also  consistent  with 
Buscheck  et  al.  (1996). 

Postulated  (based  on  the 
standard  assumption  that  a*  ~ 
10%  plume  length). 

Postulated  (based  on  the 
standard  assumption  that  oty  ~ 
10%  x  ax). 

Postulated  (assuming  source 
became  active  in  late  1960s). 


IT  (1996). 
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these  high  permeability  zones  were  found  to  be  only  on  the  order  of  2  ft  in  thickness,  so 
thecontinuity  of  such  features  over  large  lateral  distances  is  questionable.  It  is  likely,  in  an  alluvial 
environment  such  as  that  of  GAFB,  that  such  high  permeability  features  would  tend  to  appear  and 
then  pinch  out  again  frequently  across  the  aquifer.  Groundwater  flow  lines  would  tend  to 
converge  into  such  features  and  then  disperse  again  once  a  pinch  out  is  encountered.  As  a  result, 
megascopically-averaged  flow  velocities  over  large  distances  in  the  subsurface  would  be  likely  to 
approach  those  suggested  by  aquifer  tests. 

Nevertheless,  to  address  this  concern,  separate  forecast  scenarios  were  run  for  the  distribution 
of  hydraulic  conductivity  given  by  aquifer  testing  (30  ft/day  to  100  ft/day),  and  for  higher  values 
representative  of  preferential  subsurface  flow  pathways  (300  ft/day  to  1,000  ft/day,  also  assuming 
a  lognormal  distribution).  The  results  are  shown  on  Figure  C-l  for  1995  and  2022.  Two 
observations  emerge  from  this  analysis:  1)  lower  hydraulic  conductivity  values  produce  more 
realistic  estimates  of  present  benzene  plume  length  (on  the  order  of  4,000  ft  for  the  50th  percentile; 
within  an  order-of-magnitude  of  the  approximate  current  downgradient  extent  of  the  GAFB  plume) 
than  do  the  higher  estimates  (on  the  order  of  30,000  ft),  and  2),  for  both  scenarios,  forecast  plume 
lengths  do  not  change  appreciably  with  time  (i.e.,  plume  has  become  stable  regardless  of  flow 
velocity).  These  results,  taken  together,  suggest  that  the  higher  flow  velocities  measured  by 
Wilson  et  al.  (1995)  do  not  necessarily  imply  that  the  GAFB  plume  will  continue  to  grow, 
assuming  an  adequate  supply  of  sulfate  and  other  electron  acceptors. 

Sensitivity  analyses  of  the  simulation  results  suggest  that  uncertainty  in  the  degradation  rate  has 
the  greatest  impact  in  plume  length,  followed  by  hydraulic  conductivity,  gradient,  dispersivity,  and 
fractional  organic  content. 

C-2.1.  Mass  Balance  Analysis  of  Plume  Stability 

As  described  above,  a  simple  plume  model  assuming  a  continuous  source  and  a  first-order 
decay  relationship  can  be  used  to  gain  insights  into  the  stability  of  a  contaminant  plume.  As 
second,  independent  means  for  analyzing  plume  stability  consists  of  a  mass  balance  analysis, 
which  attempts  to  quantify  the  rates  of  contaminant  input  into  groundwater  (in  this  case,  through 
NAPL  dissolution),  and  the  rate  of  contaminant  loss  (through  degradation). 

The  rate  of  contaminant  (i.e.,  dissolved  benzene)  input  into  the  groundwater  can  be  estimated 
by  integrating  dissolution  fluxes  calculated  in  Appendix  B  (i.e.,  Equation  B-3)  across  the  estimated 
surface  area  of  the  NAPL  pool.  Assuming  a  solute  flux  of  1.82  x  10~3  g/m2/day  (the  median  value 
yielded  in  the  Monte  Carlo  dissolution  analysis)  and  a  contact  area  of  1 ,860,000  ft2,  the  total  mass 
of  benzene  dissolving  into  groundwater  would  be  equal  to  approximately  310  g/day.  It  is  highly 
likely  that  this  number  is  extremely  conservative  because  1)  dissolution  fluxes  out  of  the  NAPL 
pool  will  decrease  over  time  and  2)  volatilization  will  deplete  benzene  and  other  volatiles  even  more 
rapidly  than  dissolution  (see  Appendix  B).  Nevertheless,  using  this  conservative  solute  flux 
estimate,  and  a  first-order  decay  coefficient  of  approximately  0.80%  day-1  (Wilson  et  al.,  1995),  a 
total  mass  of  39.2  kg  of  dissolved  benzene  in  groundwater  would  be  needed  to  maintain  a  net  mass 
balance  (assuming  a  uniform  degradation  rate). 

Assuming  a  simplified  uniform  square  plume,  mean  concentrations  of  dissolved  benzene 
associated  with  this  mass  are  shown  on  Table  C-2  over  a  variety  of  spatial  scales.  Variability 
includes  different  plume  depths  and  area  size  factors  relative  to  the  NAPL  pool  area  (assuming  a 
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uniform  porosity  of  25%).  In  general,  given  the  observed  spatial  extent  and  depth  of  the  benzene 
plume  at  GAFB,  and  the  observed  benzene  concentration  distributions,  the  results  suggest  that 
assertion  of  plume  stability  appears  to  be  quite  reasonable  purely  from  a  mass  balance  viewpoint. 


Table  C-2.  Mean  concentrations  of  BTEX  implied  by  mass  balance  arguments. 


Plume  depth 
(m) 

Size 

factor 

Contaminated 

volume 

(m3) 

Implied  mean 
concentration 
(g/m3) 

Implied  mean 
concentration 
(ppb) 

0.1 

1 

4.32  x  103 

9.08 

9,080 

0.5 

i 

2.16  x  104 

1.82 

1,816 

1 

i 

4.32  xlO4 

9.08  x  10'1 

908 

2 

i 

8.64  xlO4 

4.54  x  10-1 

454 

5 

1 

2.16  x  10s 

1.82  x  10'1 

182 

0.1 

2 

8.64  x  103 

4.54 

4,540 

0.5 

2 

4.32  x  104 

9.08  X 10'1 

908 

1 

2 

8.64  xlO4 

4.54  xlO'1 

454 

2 

2 

1.73  xlO5 

2.27  x  lO'1 

227 

5 

2 

4.32  x  105 

9.08  x  10'2 

91 

0.1 

4 

1.73  x  104 

2.27 

2,270 

0.5 

4 

8.64  xlO4 

4.54  x  10'1 

454 

1 

4 

1.73  x  10s 

2.27  xlO’1 

227 

2 

4 

3.46  x  10s 

1.13  x  101 

113 

5 

4 

8.64  x  10s 

4.54  x  lO’2 

45 

C-3.  Reactive  Transport  Analysis  of  Geochemical  Indicators 

Geochemical  indicators  of  bioattenuation  have  been  measured  in  groundwater  at  GAFB  as  part 
of  the  natural  attenuation  study  (IT,  1996).  These  parameters  include  dissolved  oxygen,  nitrate, 
sulfate,  manganese,  iron,  methane,  bicarbonate  (as  total  alkalinity),  pH,  and  Eh.  Different 
biologically-mediated  oxidation  reactions  may  be  associated  with  each  of  the  principle  electron 
acceptors  available  (i.e.,  oxygen,  nitrate,  sulfate,  manganese,  iron).  For  example,  Bradley  et  al. 
(1992)  found  clear  evidence  of  biodegradation  of  JP-4  components  through  denitrification 
processes  in  a  shallow  aquifer.  To  gain  insights  into  the  integration  of  various  biogeochemical 
regimes  with  transport  processes  in  groundwater  at  GAFB,  a  semi-analytical  reactive  transport 
model  was  employed  to  simulate  the  chemical  evolution  of  the  impacted  aquifer  over  time.  The 
model  used  for  this  purpose  couples  an  analytical  solution  to  the  advective-dispersive  transport 
equation  with  a  geochemical  speciation  model  (McNab,  1996).  The  geochemical  speciation  model 
uses  mass  conservation  and  thermodynamic  constraints  to  estimate  changes  in  water  chemistry 
parameters  resulting  from  equilibration  reactions.  The  analytical  transport  model  uses  the  two- 
dimensional  solution  of  Domenico  (1987)  to  simulate  the  movement  of  a  contaminant  undergoing 
decay  with  a  continuous  line  source  approximation. 
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Using  a  set  idealized  hydrologic  and  chemical  parameters  characteristic  of  the  GAFB  site 
(Table  C-2),  the  evolution  of  the  local  geochemical  environment  in  association  with  the  degrading 
hydrocarbon  plume  was  simulated  over  a  period  of  10,000  days  (approximately  27  years);  a  period 
of  time  roughly  consistent  with  the  release  history  at  the  site.  For  simplicity,  the  dissolved 
hydrocarbon  plume  was  represented  chemically  by  toluene,  C7Hg,  with  a  constant  source  area 
concentration  of  40  mg/L  (depth-averaged). 

Comparisons  between  measured  indicator  data  and  simulated  plume  evolution  are  shown  on 
Figures  C-2  through  C-l  1  for  the  various  parameters  (groundwater  flow  direction  is  approximately 
north).  In  general,  model  predictions  are  qualitatively  consistent  with  observed  changes  in 
indicator  parameters  in  association  with  the  hydrocarbon  plume  (Fig.  C-2).  These  changes  include 
locally  depleted  dissolved  oxygen,  nitrate,  and  sulfate  concentrations  (Figs.  C-3  through  C-5, 
respectively),  local  generation  of  methane  (Fig.  C-6),  local  mobilization  of  manganese  and  iron  due 
to  biodegradation/thermodynamic  undersaturation  of  manganese  and  ferric  iron  oxides  (Figs.  C-7 
and  C-8,  respectively),  locally-depleted  Eh  (Fig.  C-9),  locally-increased  inorganic  carbon  as 
bicarbonate  as  a  result  of  hydrocarbon  mineralization  (Fig.  C-10),  and  a  local  decrease  in  pH 
(Fig.  C-l  1).  These  simulation  results  suggest  in  a  qualitative  manner  how  observed  spatial 
changes  in  inorganic  geochemistry  noted  at  GAFB  may  be  explained  in  a  consistent  manner  by 
bioattenuation  reactions. 


Table  C-3.  Reactive  transport  modeling  parameters. 


Pore  velocity 

2.2  x  lO-2  m/day 

Longitudinal  dispersivity  (postulated)1 

80m 

Transverse  dispersivity  (postulated) 

10  m 

Source  concentration 

40  mg/L 

Source  width 

50m 

Mean  toluene  half-life  (postulated)2 

300  days 

Background  pH 

7.6 

Background  Eh 

770  mV 

Background  02 

8  mg/L 

Background  HC03‘ 

230  mg/L 

Background  N03" 

2  mg/L 

Background  S042’ 

170  mg/L 

Differences  in  major  hydrocarbon  geochemical  indicator  parameters  between  the  hydrocarbon 
plume  and  background  groundwater  composition  are  shown  on  Figure  C-l 2  for  both  measured 
data  and  simulation  results.  For  both  the  field  data  and  the  model  results,  the  hydrocarbon  plume 
is  defined  for  this  analysis  as  characterized  by  water  composition  with  total  BTEX  greater  than 

1  ppm,  while  background  is  defined  as  total  BTEX  less  than  0. 1  ppm. 

'Large  dispersivity  values,  in  comparison  to  the  scale  of  observation,  were  chosen  for  this  analysis  as  conservative 
assumptions  to  account  for  the  presence  of  preferential  flow  pathways. 

2 

Conservative  assumption. 
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Several  important  observations  emerge  from  these  data:  1)  approximately  80  percent  of  the 
observed  depletion  of  BTEX  between  plume  conditions  and  background  conditions  can  be 
accounted  for  stoichiometrically  in  the  accumulation  of  inorganic  carbon  (bicarbonate),  which  is 
indicative  of  hydrocarbon  mineralization,  2)  observed  net  changes  in  average  concentration  of  the 
major  electron  acceptors  accounts  for  approximately  2.7  times  the  observed  generation  of 
bicarbonate.  This  indicates  that  biologically-mediated  oxidation  of  the  hydrocarbons  (primarily  by 
sulfate  reduction)  can  easily  account  for  the  observed  transformation  of  the  hydrocarbons  into 
inorganic  carbon  (the  apparent  loss  of  a  significant  portion  of  the  inorganic  carbon  may  be  due  to 
offgassing  of  C02),  and  3)  when  comparing  observed  changes  in  indicator  parameters  to  those 
predicted  by  the  model,  it  is  apparent  that  the  model  is  actually  conservative  in  estimating  the 
amount  of  dissolved  hydrocarbon  that  is  biodegraded  in  the  center  of  the  plume  (e.g., 
methanogenesis  and  sulfate  reduction  are  underestimated;  inorganic  carbon  production  is 
underestimated). 
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Figure  C-1.  Forecast  idealized  benzene  plume  lengths  for  low  and  high  hydraulic  conductivity  test  cases. 
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Figure  C-2.  Dissolved  hydrocarbon  -  measured  total  BTEX  (left,  mg/L)  and  simulated  toluene  (right,  mol/L). 


-  measured  (left,  mg/L)  and  simulated  (right,  mol/L). 


Figure  C-4.  Nitrate  •  measured  (left,  mg/L)  and  simulated  (right,  mol/L). 


Figure  C-5.  Sulfate  -  measured  (left,  mg/L)  and  simulated  (right,  mol/L). 


e  C-6.  Methane  -  measured  (left,  micrograms/L)  and  simulated  (right,  mol/L). 


Ferrihydrite  dissolution  favored. 


Figure  C-8.  Iron  -  measured  (left,  micrograms/L  in  solution)  and  simulated  thermodynamic  saturation  with  ferrihydrite  (right). 


red  (left)  and  simulated  (right),  millivolts. 


e  C-10.  Bicarbonate  -  measured  (left,  mg/L  CaC03  bicarbonate  alkalinity)  and  simulated  (right,  mol/L). 
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Appendix  D 


Inhalation  Risk 

The  inhalation  risk  at  the  surface  was  evaluated  using  a  one-dimensional  diffusive  transport 
model  for  flux  through  the  vadose  zone.  The  volatile  flux  of  chemicals  in  a  1-D  porous  media  have 
been  assessed  by  Jury  et  al.  (1983),  and  Jury  et  al.  (1990)  using  an  analytical  modeling  approach. 
Factors  taken  into  account  included  degradation,  as  well  as  an  effective  water  solute  velocity,  the 
latter  which  may  be  used  to  assess  sites  with  either  leaching,  evapotranspiration,  or  no  water 
movement  through  the  subsurface.  At  GAFB,  evapotranspiration  is  expected  to  far  exceed  that  of 
the  yearly  recharge.  However,  due  to  the  extreme  depth  of  the  groundwater  zone,  at  38  meters 
below  the  surface,  evapotranspiration  is  not  a  likely  factor  in  the  transport  of  volatile  components. 
Therefore,  the  version  of  the  Jury  model  presented  here  does  not  take  this  process  into  account. 

According  to  Jury  et  al.  (1990),  the  cumulative  flux  at  the  surface,  assuming  a  zero  boundary 
layer  thickness,  can  be  expressed  as: 


Kb(°°)  =  C0H 
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where  Vcb(°°)  is  the  cumulative  vapor  flux  (g-cm),  Ca  is  the  concentration  of  the  volatiles  diffusing 
from  the  separate  phase  FHC  product  (g/L),  L  the  depth  to  the  NAPL,  W  the  thickness  of  the 
NAPL,  A  the  degradation  rate  (1/sec),  and  DE  is  the  effective  diffusion  coefficient  as  defined  by: 

+  G]mD, 

°E  ~  «!(p»/^+e,+e^)  <D'2) 

where  the  equation  parameters  are  as  defined  in  Appendix  B.  D,  is  the  diffusion  coefficient  of 
volatile  component  in  water. 

The  target  risk  level  associated  with  the  above,  is 
Cg  ■  CPF  IR  EF-  ED 

TR  ~  BW~AT^365dfJr  (D‘3) 

where  CPF  refers  to  the  cancer  potency  factor,  IR  is  the  inhalation  rate  of  air,  EF  is  the  exposure 
frequency,  ED  is  the  exposure  duration,  BW  is  adult  average  body  weight,  and  AT  is  the  averaging 
time,  as  determined  in  Dooher  et  al.  (1996). 
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Since  Equation  D-l  reflects  the  cumulative  flux  over  an  infinite  time  period,  it  must  be  divided 
by  the  averaging  time,  AT,  as  well  as  a  mixing  zone  height,  d,  taken  here  to  be  2  m,  to  give  the 
yearly  average  concentration.  This  gives: 


V‘H 

AT-d 


(D-4) 


A  Monte  Carlo  simulation,  utilizing  Crystal  Ball  in  a  manner  similar  to  that  described  for  plume 
length  simulation  in  Appendix  C,  and  addressing  all  the  above  aspects,  produced  two  very 
different  results  based  on  differing  degradation  rates.  The  median  risk  level  was  found  to  be  a 
4  x  1CT5  excess  lifetime  cancer  risk,  with  the  90th  percentile  equal  to  9.8  x  10~4,  using  a  first 
order  degradation  rate  equivalent  to  that  used  in  the  groundwater  plume  stability  analysis 
(9.95  x  10“8/sec). 

As  degradation  rates  in  the  vadose  zone  are  not  considered  to  be  oxygen  limited,  due  to  the 
high  potential  flux  from  the  atmosphere,  and  since  the  analysis  deals  with  volatiles  sorbing  onto  the 
surface  film  surrounding  soil  particulates,  oxygen  is  much  more  freely  available  than  found  in 
groundwater.  Therefore,  the  rate  mean  and  standard  deviation  as  relevant  to  groundwater,  were 
both  increased  by  an  order  of  magnitude  in  order  to  simulate  the  lower  half-lives  found  in  a  more 
aerobic  environment  than  groundwater.  Proceeding  as  before,  this  produces  a  median  risk  level 
well  below  10”6,  with  a  90th  percentile  level  of  4.5  x  1 0  \  The  largest  contribution  to  variance  in 
the  above  analysis  is  the  degradation  rate,  followed  by  water  content  of  the  soil.  Both  cumulative 
distributions  are  shown  on  Figure  D-l.  The  distributions  not  listed  elsewhere  in  these  appendices 
are  found  in  T able  D- 1 . 

It  should  be  noted  that  the  above  calculations  are  for  a  1-D  flux  under  the  influence  of  a  zero 
concentration  boundary.  In  reality,  the  risk  is  much  lower  due  to  the  near  zero  flux  condition 
imposed  by  the  asphalt/concrete  cap.  Hence,  these  results  should  be  considered  very  conservative. 


Table  D-l.  Inhalation  risk  modeling  parameters. 


Parameter 

Probability  distribution 

Basis 

Body  weight  (BW) 

Normal.  Mean  =  71  kg;  standard 
deviation  =  15.9  kg 

Dooher,  1996 

Exposure  duration-Westem  U.S. 
(ED) 

Lognormal.  Mean  =  3.49  years; 
standard  deviation  =  6.84  years 

Dooher,  1996 

Exposure  frequency  (EF) 

Triangular.  180,340,365 

Dooher,  1996 

Inhalation  rate  (IR) 

Lognormal.  Mean  =  12.06  L/day; 
standard  deviation  =  2.55  L/day 

Dooher,  1996 

Cancer  potency  factor  (CPF) 

Point.  0.10  kg  day/mg 

Dooher,  1996 

Averaging  time  (AT) 

Point.  70  years 

EPA 

Volumetric  water  content  (0,) 

Lognormal.  Mean  =  0.04; 
standard  deviation  =  0.04 

Based  on  site  data  (IT,  1996) 

Volumetric  air  content  (0g) 

Function  of  water  content  and 
porosity 

Based  on  site  data  (IT,  1996) 
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Figure  D-1.  Cumulative  distribution  of  lifetime  excess  cancer  risk  for  benzene  inhalation  at  GAFB,  with  two  different  degradation  rates. 


